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Limit-Cycle Stability Study of a 


Feedback Control System by a New 


Describing-Function Technique 


H. J. HARRINGTON 


ASSOCIATE MEMBER AIEE 


A NEW describing-function technique 
is used in this paper to show the 
ects of several gain-controlling non- 
earities on the limit-cycle stability of an 
craft roll-stabilization system. The 
w method leads to a stability diagram 
hich is divided not into the usual two 
sions of ‘‘stable’’ and “‘unstable’’ opera- 
n but rather into numerous regions each 
which is associated with a particular 
pe of limit-cycle oscillation or singular 
int. Such a diagram gives considerable 
ight into the mechanism of oscillations 
en observed in physical equipment. 
ry satisfactory agreement with analog 
mputer studies is demonstrated for the 
w technique. 
Whe case of simultaneous limiting in 
h of two feedback loopsiseasily handled 
the system studied herein. A logical 
athesis procedure for deliberately in- 
ducing a new nonlinearity in order to 
mteract an existing one is also pro- 
paper’ treating the stability of linear 
ems furnished the impetus for this 
tk; another paper? concerned itself 
h some of the same problems which 
treated here but the approach and 
erpretation of results were considerably 
erent. 
Although the technique developed here 
y be employed for any nonlinearity for 
ch a describing function can be derived, 
appears to offer interpretative ad- 
atages over the more conventional ap- 
3 only in the case of gain-controlling 
ilinearities; that is, nonlinearities 
ich do not introduce phase shift to 
als passing through them. 
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Nomenclature 


A=peak amplitude of a sinusoidal oscilla- 
tion in roll rate 

d=value of threshold in rate gyroscope 

Em=Ma—M;, error in normalized rolling 
moment 

f=viscous damping coefficient of air frame 

I=moment of inertia of air frame about 
roll axis 

t=control surface incidence 

t¢=control surface incidence command 

J=V—-1 

Ka=Mi/i, aerodynamic effectiveness of 
control surface 

Kp=position feedback gain 

K,=rate feedback gain 

Ly=position-gyro limit 

L,=rate-gyro limit 

Mg=rolling moment (normalized to [) 
caused by external disturbance 

M,=rolling moment (normalized to I) 
‘caused by control surface incidence 

Np=equivalent gain, or describing function, 
of nonlinearity in position gyro 

N, =equivalent gain, or describing function, 
of nonlinearity in rate gyro 

s=Laplace operator 

Ti, T2=constants of control servo 

aq=d/A, amplitude parameter for threshold 
describing function 

a =L,/A, amplitude parameter for limiting 
describing function 

y=L,/L,, position-to-rate limit ratio 

¢=roll-position angle 

Q=yw, normalized frequency 

w=circular frequency 


Description of System Studied 
The system to be studied, which is 
shown in block diagram in Fig. 1, is a 


typical aircraft roll-stabilization system 
used to maintain air-frame roll angles near 


Mat 


Fig. 1. Block dia- 
gram of an_ aircraft 
roll-control system 
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zero in the presence of external disturb- 
ances. The system works as follows. A 
disturbing moment Ma produces a 
moment error which in turn causes the air- 
frame toroll. The instantaneous roll rate 
and position are measured by gyroscopes 
whose output signals are suitably scaled 
by the gains K, and K,, added, and sent 
to a control surface servo as an incidence 
command 7%. The output incidence of the 
control servo then gives rise to a control 
moment M;, which tends to reduce the 
moment error and thus stabilize the air- 
frame. 

This seemingly simple system is com- 
plicated by the presence of limiting in, or 
immediately following, the position gyro 
and by both limiting and threshold in the 
rate gyroscope. In addition, the aero- 
dynamic gain K, is subject to rather wide 
variations from nominal because of un- 
predictable changes in aircraft flight con- 
ditions; good stability margins with 
respect to the gain Kz must therefore be 
attained to prevent either divergent os- 
cillations or a hard limit cycle of cata- 
strophic proportions. 

In the study which follows, it will be 
assumed that the linearized transfer func- 
tion of both gyroscopes is unity. Al- 
though this might well be a good assump- 
tion, the restriction is not at all necessary 
and is made only for the sake of brevity 
as is the omission of compensating net- 
works in either of the two feedback paths. 
The nominal system parameters for the 
rest of the system are as follows: 


f/I=0.67 per second 
T, =0.04 second 


Ka d P 
Kp Ln? general 
T2=0.05 second Ke 1h 


(1) 
Generalized Stability Chart 


The characteristic equation F(s) of the 
system of Fig. lis 


re)=Test+( ret ret) 


(147-4 )4( xan toot 


KaKypNp (2) 


AIRFRAME 


SERVO 
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EFFECTIVE GAIN, Nr(oCy) 


EFFECTIVE RATE GAIN, KaKrNr 


40 60 80 100 
EFFECTIVE POSITION GAIN, KaXpNp 


Fig. 2. Stability diagram for the roll-control system 


For the present the nonlinearity terms 
N, and N; are assumed to be simple con- 
stants. Setting equation 2 equal to zero 
and replacing s by jw gives an equation 
which represents the marginally stable 
case if real roots in w exist: 


ties (: +E +KoK Ny 


jo| (n+-72%) ot (xaxv-+4) | =0 


(3) 


But for equation 3 to be valid, both the 
real and imaginary parts must be identi- 
cally zero; therefore, this equation may be 
replaced by two separate equations from 
which the system gains can be solved 


KyKpN> = ( ES rat 1475 ot 


(real part) (4) 


KaK;,N; = (rien Nef 
(imaginary part) (5) 


For the particular parameter values 
in equation 1, these equations become 


KaKpNp =(—0.0016w2+ 1.033)? (6) 


Fig. 3. 


K,K,N, =0.0511w? —0.67 (7) 


A plot of equations 6 and 7 with w 
treated as a real-value parameter is shown 
in Fig.2. By applying Hermite’s criteria! 
(or any other convenient method), it can 
be shown that the region within the 
closed boundary is stable while that 
outside the boundary line is unstable. 
This situation is approximately true even 
when the system is nonlinear if the effec- 
tive gains of the nonlinearities NV, and N, 
are explicitly included in the gain co- 
ordinates as shown in Fig. 2. 


Limit-Cycle Stability Chart 


For better interpretation of the non- 
linear cases, it is desirable to have the 
stability plot in terms of just the linear 
gains K,K, and K,K,, and show not only 
the wholly stable and unstable regions 
but also the various limit-cycle oscilla- 
tion conditions which might arise. To 
illustrate this, the case of a single non- 
linearity, namely, threshold in the rate 
gyro, will be examined. 


It will be assumed that the effective 


Fig. 4 (left). Limit- 
cycle stability dia- 


gram for case of 
threshold in rate 


gyroscope 


—_ 


LINEAR RATE GAIN, KaKr 


Fig. 5 (right). Phase- 


plane sketches show- 


130 Harrington—Stability Study of Feedback Control by a New Technique 


ing typical trajec- 

tories associated 

with the seven re- 
gions of Fig. 4 


OSCILLATION AMPLITUDE RATIO, OCa=/A 


Describing function for threshold 


0.6 0.8 1.0 1.2 


gain of the threshold element J, is git 
by its describing function*® during 
oscillation. A plot of N, as a function ¢ 
normalized roll-rate amplitude ag is sho 
in Fig.3. Note that if the roll-rate amy 
tude is equal to or less than the deadspe 
that ig ag>1, the gain NV, is zero. As 
roll-rate amplitude approaches infini 
however, the grain N, approaches ut 
and the system should behave linearly 
The consequences of this nonlinea 
may be seen by examining Fig. 4. H 
the solid curve is identical to that of 
2 while the broken lines, which re 
the presence of the threshold, divide tt 
over-all chart into regions of vary 1 | 
degrees of stability. In general, linear 
gain co-ordinates appearing as a pot 
Fig. 4 will be transformed into a ve 
line in Fig. 2. Since N, is by defin 
unity, this line will be located at the sama 
abscissa and will descend from the poit 
ordinate (N,=1) to the position-gain axq 


=a 


ee aS ee" 


STABLE 
LIMIT CYCLE 


>18 Lees 
<252 UNSTABLE 
: LIMIT CYCLE 


P 


PLATT s, 


Het : ph 

. \\\ 

3 \X $ 
iS 


LINEAR POSITION GAIN, KaKp, db 


If the linear-gain co-ordinates 
I in region I of Fig. 4, no amount of 
tenuation because of N,(e@a) can cause 
€ corresponding gain ordinate in Fig. 2 
) fall in the unstable region. Thus, for 
gion I, any initial disturbance must 
timately die out as is depicted by the 
brresponding phase-plane sketch of F ig. 
~ (The phase-plane sketches of ise, 
€ for purposes of qualitative description 
ly since a polydimensional phase space 
fequired to accurately describe higher- 
der systems.) 
Tf the linear-gain co-ordinates fall in 
ther regions II or III, it is possible to 
cite an unstable limit-cycle oscillation. 
is behavior may be predicted by noting 
om Fig. 2 that the system is stable when 
rest (V,=0), and-unstable when excited 
ficiently to cause the gain K,K,N, to 
ceed the stable-region ordinate. 
Region IV of Fig. 4, which is the most 
obable region of operation for the 
tem in question, leads to a soft, stable 
mi For examination of Figs. 2 
d 3 shows that with the system at rest 
e effective rate gain is zero and the sys- 
m is unstable; but, as the amplitude 
the resulting divergent oscillation 
lds up, aq decreases and N,(a,) in- 
eases until the gain reaches the lower 
ability boundary of Fig. 2. Further 
sreases in the oscillation amplitude force 
e effective gain ordinate into the stable 
zion with the result that the excess 
iplitude is damped out. Obviously the 
quency of oscillation must lie in the 
ige of 3.6<w< 18 radians per second. 
Both a soft, stable limit cycle and a 
sher-amplitude unstable oscillation are 
ssible if the system gains are adjusted 
lie in region V of Fig. 4. Repetition 
the preceding logic will show that the 
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Fig. 6 (left). Loga- 1.0 
rithmic plot of limit- 


cycle stability dia- 
gram showing over- 


lay of describing 


function for rate- Tt 
e 
gyro threshold Jas 
= 
: 
= 04 
: 
0.2 


OSCILLATION FREQUENCY, U, RADIANS /SECOND 


> 


0 20 40 


0 
80 100 120 140 160 180 200 
LINEAR POSITION GAIN, Kak, 


Fig. 7. Amplitude and frequency of limit-cycle oscillations caused by 
rate-gyro threshold for variations in position gain K.K,, while holding 
rate gain constant at K,K,=20 


Stable oscillations as determined by describing functions 


—-~ Unstable oscillations as determined by describing functions 
O Points obtained by analog computer 


stable oscillation will occur at a frequency 
corresponding to the lower stability 
boundary of Fig. 2 while the unstable 
limit cycle will be at one of the higher 


frequencies of the upper stability bound-— 


ary. 

Finally, regions VI and VII represent 
totally unstable conditions for the case of 
rate-gyro threshold just as they do for a 
linear system since for these gains no 
realizable value of N;(ag) can lead to the 
stable region of Fig. 2. As shown in Fig. 
5, the phase-plane sketch characterizing 
these regions consists of typical tra- 
jectories near an unstable singular point. 

The preceding discussion has shown 
that the stability plot for the linear system 
(Fig. 2, N,=N,=1) consists of only two 
separated regions which define the sys- 
tem as being either stable or unstable. 
The corresponding limit-cycle stability 
chart for the case of rate-gyro threshold, 
however, is divided into seven separated 
regions which represent five distinct types 
or degrees of stability. 


Quantitative Limit-Cycle Data for 
Threshold 


So far, only the nature of possible limit- 
cycle oscillations caused by threshold have 
been discussed. Now the amplitude and 
frequency of the oscillations will be 
quantitatively determined and compared 
with the results of analog computer tests. 

The normalized amplitude ag of the 
limit cycle can always be determined by 
computing the attenuation necessary for 
oscillations from Fig. 4 and then referring 
back to Fig. 3 to see what ag value gives 
this attenuation. A more convenient and 
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elucidating method, however, is to replot 
the portions of the stability chart and 
describing function of interest to a log- 
arithmic scale; the describing-function 
graph can then be used as an overlay and 
the gain division indicated previously be- 
comes simple graphical subtraction. 

As an illustration, portions of Figs. 3 
and 4 are replotted in terms of decibels 
(db) in Fig. 6. Here, the describing- 
function graph is drawn as N,(ag) versus 
N, with ag as a parameter; since the posi- 
tion gyro is presently assumed linear, NV Fs 
is identically unity (zero db) for all values 
of ag. Two separate overlays of the 
describing function are shown on the sta- 
bility chart of Fig. 6 as they would appear 
for the following gain conditions 
Case A: K,K,=26 db; K,K,=35.6 db 
Case B: K,K,=31 db; K,K,=43 db 
For each of these system gains, the oscilla- 
tion amplitude and frequency may be 
immediately determined by noting the 
values of ag and w corresponding to the 
intersections of the two curves. Note 
that for case B, both a stable and un- 
stable solution exists for respectively the 
lower and upper intersection. 

The accuracy with which the amplitfide 
and frequency of limit-cycle oscillations 
can be predicted is illustrated in Fig. 7 
which compares the results of analog com- 
puter studies and the proposed analytical 
technique. The analog setup employed a 
conventional biased-diode circuit for the 
threshold simulation; all harmonics 
from this circuit were allowed to circulate 
freely around the loop. The system rate 
gain was held at KzK,=20 (26 db) and the 
position gain then varied to obtain the 
amplitude and frequency points and 
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10 Fig. 8 (left). De- 
scribing function for 
limiting 

08 
3 Fig. 9 (right). De- 
ae scribing - function 
e overlay for simulta- 
neous limiting in 
position and rate 
f gyros. Each curve 


is associated with a 


0 0.2 0.4 0.6 0.8 


OSCILLATION AMPLITUDE, oC, =L:/A 


curves shown in Fig. 7. The darkened 
triangles appearing on the analytically 
determined curves correspond to the 
similarly marked point in Fig. 6, case A. 

It may be noted from Figs. 4 and 5 that 
the gain conditions K,K,=20, 156<KiKy 
<167 corresponds to region V wherein 
both stable and unstable oscillations are 
predicted. Both types were verified by 
the analog computer study and, as shown 
in Fig. 7, the proposed analytical method 
is capable of determining the amplitude 
and frequency of either type oscillation 
with quite suitable engineering accuracy. 
(It would appear, however, that the same 
reservations concerning higher harmonics 
applicable to conventional describing- 
function theory? are equally applicable to 
the new method.) 


Rate-Gyro Limiting 


The effects of limiting in the rate gyro 
can be determined by employing the same 
procedure just discussed for threshold 
after substituting the now-appropriate 
describing function for limiting shown in 
Fig. 8. Such an investigation revealsthat 
the limit-cycle stability diagram re- 
mains divided into the same seven regions 
shown in Fig. 4. In all but region V, 
however, the stable and unstable oscilla- 
tions shown in Fig. 5 become inter- 
changed in going from threshold to limit- 
itt; in region V the frequencies of the 
two limit cycles become interchanged. 

Analog-computer studies of the limit 
cycles resulting from rate-gyro limiting 
were also conducted. The analytical re- 
sults were verified by this means to within 
essentially the same accuracy as that 
shown in Fig. 7 for threshold. 

Investigation of the rate-gyro limiting 
case leads to a very interesting conclusion. 
The whole of region IV, Fig. 4, is plagued 
by a hard limit-cycle oscillation and yet 
this is the only region of practical opera- 
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single value of nor- 
malized frequency 
Q= yw 


1.0 1.2 14 


tion from a linear standpoint. Since the 
threshold of a physical rate gyro usually 
increases in direct proportion to the maxi- 
mum range, or limit value, it is usually 
not wise to specify a very large limiting 
value in an attempt to rationalize the 
limit cycle out of existence. What, 
then, can be done to circumvent this 
problem area? One answer is to delib- 
erately introduce another nonlinearity 
into the system. 


Rate- and Position-Gyro Limiting 


It is obvious that position-gyro limiting 
would cause the describing-function loci 
of Fig. 6 to swing towards the left, for 
the value of N, would then become less 
than unity. This suggests that operation 
in region IV, Fig. 4, might be possible 


without the potential hard limit cycle if. 


position-gyro limiting were employed. 

To investigate this possibility, assume 
that an oscillation in roll exists such 
that 


$=A sin ot (8) 


s= [dan—< cos wt (9) 


If now the rate signal is limited to the 
value +L,, and the position signal is 
limited to +L,, the describing function 
N(a) for each limiter will be as shown in 
Fig. 8; here, a is defined for the rate 
limiter as 


ate 10 
ay) =1 ( ) 
and for the position limiter, 

_Lyp LLy 
ave 7, Caerernge (11) 


From equations 10 and 11, a new over- 
lay chart of N,(ai) versus Ny(a2) may 
be constructed with loci for several values 
of normalized frequency, 2= yw, as shown 
in Fig. 9. This overlay can be used in 
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RATE LIMITER GAIN, Ne(cC), db 


eal 
POSITION LIMITER GAIN, Np( OC; 2 ), db 


conjunction with a linearized-gain § 
bility plot just as was done in Fig. 6 
cept now only those intersections o¢ 1 
ring at like frequencies are valid, T 
region can be found in which no 
intersections occur, then no limit ecyel 
due to gyro limiting are possible in thi 
region. J 
The determination of such a regi 
demonstrated in Fig. 10 which sh 


‘few values of the N;(a1) versus Nj(e 


overlay drawn on the stability chart? 
the special case of y=L,/L,;=1/5 secon 


32 


LINEAR RATE GAIN, KaKr, db 


4 
28 
LINEAR POSITION GAIN, KaKp, db 


Fig. 10. Logarithmic plot of stability di 
and describing-function overlay for 
taneous gyro limiting in the ratio of 7 


—— — Separates regions of stable and hare 
cycle operation 

A, Boundary solutions obtained by di 

ing Function 

LJ Solutions obtained from analog co 


28 32 36 40 44 48 
LINEAR POSITION GAIN, KaKp, db 


11. Logarithmic plot of stability diagram 
showing synthesis procedure 


e that if the describing-function over- 
were moved up from the gain value 
wn, no more valid (like-frequency) 
srsections would occur. But if the 
a were lowered, two proper intersec- 
as would result: one corresponding to 
ard, unstable limit cycle of frequency 
10.6 radians per second, and the other 
able oscillation of somewhat greater 
plitude and lower frequency. It is 
y apparent that as the system gain 
roaches a value for which valid inter- 
ions are only marginally possible, the 
limit cycles coalesce to destroy one 
ther. Such a gain, therefore, repre- 
is one point on a boundary line which 
arates the stability chart into regions 
fable and potentially unstable (hard 
t cycle) operation. 

he heavy dashed line of Fig. 10 shows 
analytically determined locus of sys- 
gains resulting in neutral limit-cycle 
ility for y=1/5 second. Also shown 
comparison are the amplitude and 
uency of several of these marginal 
ations as determined both ana- 
ally and by the analog computer. 
= the excellent agreement in all 
ects. 


mthesis Procedure 


aving developed the tools necessary 
andle the nonlinearities in the system, 
ose now that it is desired to determine 
best gains and gyro limit values to 
fy the following design requirements: 


The feedback gain ratio shall be ad- 
sd such that K,/Kp=0.15 second. 


The linear-system gain margin with 
ect to K, shall be 10 db. 


3. The hard limit-cycle gain margin with 
respect to Kg shall be 6 db. 


4. The system shall be capable of measur- 
ing roll angles up to +25 degrees. 


5. _The amplitude of soft, limit-cycle 
oscillations due to rate-gyro threshold shall 
be held to a minimum. 


The first requirement is based on a 
linear root locus study in which the 
amount of rate-to-position feedback re- 
quired to provide an adequate transient 
response in the linear mode of operation 
was determined. 

In the second and third requirements, 

the gain margin is expressed in terms of 
Ka since, as mentioned previously, this 
gain is subject to rather wide and un- 
predictable variations. The limit-cycle 
gain margin is specified somewhat lower 
than that for the linear case since not only 
the system gain but also the level of 
excitation existing simultaneously must 
be favorable before oscillations can de- 
velop. The relative weighting of the two 
gain margins is largely a matter of judge- 
ment. 
The fourth and fifth requirements state, 
in essence, that the position-gyro limit 
be set at 25 degrees. For if this limit 
were higher than necessary, the rate-gyro 
limiter would also have to be higher to 
provide a given limit-ratio value y, and 
this in turn would increase the threshold 
level of the rate gyro because of physical 
construction considerations. 

The synthesis procedure is shown in 
Fig. 11. First, the linear frequency locus 
is plotted as curve A. Next, curve B 
is established by projecting each point of 
curve A back and down 10 db; obviously 
curve B represents a boundary outside of 
which the linearized gain margin with 
respect to Kg is less than 10 db. The 
condition for K,;/K,=0.15 is depicted 
as the straight line C. Thus the point 
X satisfies the first two requirements 
previously listed and represents the 
desired gain condition. 

By virtue of requirement 3, the point Y 
may be determined as falling on the limit- 
cycle stability boundary; therefore the 
describing-function overlay, curve D, 
must furnish a marginally stable solution 
when its zero-zero db point is aligned with 
Y. This marginally stable solution is 
shown as point Z in Fig. 11 and occurs for 
w=9.8 radians per second and Q=yw= 
3.4. From these data, the limit ratio y 
may be solved for as 

Jie RG 


3.4 
Ve Deo second (12) 


Now from requirement 4 and equation 12, 
the desired value of rate limit, L,, be- 
comes 


Ly (25X9.8 
L,=— = ENR) =72 degrees per second 


Y 3.4 
(13) 


In summary, the system which meets 
the specified requirements has the follow- 
ing parameters: 


K,K,=16.6 db=6.8 per second 
K,K,=33.2 db=46 per second squared 
Ly=25 degrees 

L,=72 degrees per second 


Conclusions 


The new describing-function technique 
presented herein appears to offer quite 
substantial advantages in interpretation 
and flexibility over the more conventional 
approach to limit-cycle analysis when the 
system nonlinearities are of the gain- 
controlling type. A slight modification? 
would allow phase-shifting nonlinearities 
to be handled also but with dubious bene- 
fits. 

The additional utility of the new 
method stems from having two generalized 
system parameters with which to work 
rather than the single gain associated 
with open-loop diagrams.* The price 
paid for this convenience is the some- 
what laborious derivation and computa- 
tion of the characteristic equation of the 
closed-loop system, particularly if the 
order of the equation is high. To obtain 
the same amount of information and in- 
sight by the conventional method would, 
however, lead to vastly more labor since 
families of open-loop frequency loci would 
be necessary to account for the second 
parameter. 

For the system analyzed in this paper, 
it was convenient to choose two gains for 
the co-ordinates of the stability diagram. 
In other problems it might be more appro- 
priate to choose one gain and some other 
system parameter, a time constant for 
example. In other words, the method is 
considerably more flexible than is demon- 
strated by the particular analysis carried 
out herein. 
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Calvin Johnson. AIEE Transactions, vol. 71, pt. 
II (Applications and Industry), July 1952, pp. 
169-81. 


4, Non-Linear Mecuanics (book), N. Minorsky. 
J. W. Edwards, Inc., Ann Arbor, Mich., 1947, pp. 
68-73. 
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Discussion 


Toshio Numakura (Hitachi Central Re- 
search Laboratory, Kokubunji, Tokyo, 
Japan): The author has solved a com- 
plicated system which contains several 
nonlinearities by using the stability dia- 
grams introduced in references 1 and 2 
of the paper. 

The new ideas added to the describing- 
function method and the works cited are 
very interesting. 

Some problems may remain unsolved 
in the last section of the paper, ‘‘A Synthesis 
Procedure,’”’ in which the author intends 
to predict the transient response of the 
nonlinear system by use of the stability 
diagram. Fig. 11 does not provide any 
indication of the phase margin or slope of 
the gain characteristic as does Bode’s 
diagram, so it is not possible to predict the 


transient response exactly. Though the 
discusser and his colleagues intended to 
find a suitable method by which the 
transient response could be predicted, 
they failed to find any not only in the case 
of the nonlinear system but for the linear 
system as well. Therefore the optimum 
parameter conditions should be found by 
using the stability diagram in conjunction 
with an analog computer. 


H. J. Harrington: The author appreciates 
the comments of Dr. Numakura and the 
resulting opportunity to elaborate on the 
paper. 

Actually, there is no intent in the paper 
to predict the transient response by means 
of the stability diagram. Thus, in the 
section referred to by the discusser, it is 
pointed out that a complementary root- 
locus analysis furnished the gain ratio value 


of K,/Kp=0.15 required to give an adeqi 
transient response. It is true tha 
linear-system gain margin of 10 db 
respect to K, is marked off on the stah 
diagram (Fig. 11), but this is mostly 
matter of convenience for presente 
since the placement of point X was 
termined by the afore-mentioned 
locus study. 7 
The difficulty experienced by Dr. N 
kura and his colleagues in determinin, 
phase margin and transient response - r 
the stability diagram is shared by 
author. This limitation is so severe 
when used alone the stability diagra 
almost useless as a linear synthesis 
As expanded by this paper, however, ¢ 
method is extremely valuable in establi 
the nature of steady-state limit cycles 
mildly nonlinear systems in which 
linear-range transient response has bi 
determined by other means. ) 


New Methods of Simplifying Boolean 


Functions 


R. L. HOWARD 


ASSOCIATE MEMBER AIEE 


CONOMY IS ONE of the most im- 
portant considerations in the design 
of any system. Maximum savings can be 
achieved in the design of switching cir- 
cuits by performing the desired logic 
in the simplest minimum form applicable 
to the equipment utilized. Boolean 
algebra is often used to express the func- 
tions of switching circuits and these 
Boolean expressions generally appear in 
the final form of a sum of products or a 
product of sums. The form that pre- 
sents the most economical application is 
determined by the equipment. It is 
therefore necessary to be able to trans- 
form and simplify Boolean functions so 
that they will appear in the predetermined 
form that will yield the maximum savings. 
Methods for performing the simplifica- 
tion and transformation of Boolean func- 
tions have been presented in previous 
publications by Karnaugh,! Quine,? and 
McCluskey. Each of these processes 
becomes extremely complicated for com- 
plex functions. Furthermore, some dup- 
lications are encountered which repeat 
information during the process that 
initially had been stated in the original 
function. 

The following discussion describes a new 
procedure of simplifying and transform- 
ing any Boolean function without some of 
the restrictions encountered in previously 
presented procedures, The original 
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Boolean function may consist of any 
number of variables and may be expressed 
in any form (sum of product, product of 
sums or a combination of both). With- 
out changing the form of the original func- 
tion the same basic rules apply when 
either simplification and transformation 
is made; only the recording of the original 
function is changed. This provides a 
systematic process for the reduction of 
static or conventional control systems and 
also an easily programmed process for 
digital computer application. 


Discussion 


For clarity the explanation of the 
method includes a sample Boolean func- 
tion. Wherever the sample is concerned 
in the explanation that portion will be 
labelled “‘sample.”” Also included wher- 
ever possible, are examples which are 
independent of the sample function. 
Each of these will be labelled ‘“‘Example.” 

The sample Boolean function is 
X=A-BICD-E+A-C-E(C+D)| __ 

; =AB|CDE+ACE(C+D)] 
where 


+ =or function 
* =AND function 


The “variables” of a function are the 
number of different letters concerned. 
The ‘“‘literals’” of a function are the in- 
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contained in the original expression 


dividual entries in the function. 
the sample function contains five 
ables and ten literals. Fig. 1 sh 
logic diagrams of the sample functie 


THE BASIC PROCEDURE OF THE MBi 


Preparation of the Function for 
Transformation and Simplification 


1. Determine the number of vari 


Sample: There are five variables conte 
in the expression (A, B, C, D. and £) 


2. Assign a digit to each variable. — 


Digits 
Sample: E D CBA ‘ 


3. The digital representation of 
product is now written as follows: 


If the literal indicates the presence « 
variable, a “‘1’’ will be placed in the 
responding digit for that term. j 
If the literal indicates the absence (negé 


of the variable, a ‘‘0” will be placed i 
corresponding digit for that term. 


If the literal does not appear in the t 
dash (—) will be placed in the correspot 
digit for that term. 


Sample: AB[CDE+ACE(C+D)] 
ED. C°B A -o8 Dale 
— ie Ot 
EDCBA  F Dag 
0 - Oe 1 eae 


4. For the conversion of the 01 


Paper 59-91, recommended by the AIEE In 
Control Committee and approved by the 
Technical Operations Department for presi 
at the AIEE Winter General Meeting, 
N. Y., February 1-6, 1959. Manus 
mitted November 3, 1958; made av: 
printing December 16, 1958. 


R. L. Howarp is with Westinghouse 
Corporation, Pittsburgh, Pa. ¢ 


ession into the sum of products 
fich will be referred to as the AND 
n): (a). Wherever an AND function is 
ated in the original expression the two 
ps it concerns will be combined. (6), 
prever an OR function is indicated in 
Joriginal expression the two terms it 
terns will not be combined. 


For the conversion of the original ex- 
sion into the product of sums (which 
be referred to as the or form): (a), 
rever an OR function is indicated in 
original expression the two terms it 
erns will be combined. (b). Wherever 
D function is indicated in the original 
ession the two terms it concerns 
not be combined. 


Symbols used for indicated combina- 

and noncombinations are a long 
(—) between two terms to indicate 
they are to be combined; and a long 
ed dash (-+-) used between two 
to indicate that they are not to be 
bined. These symbols will vary in 
to correspond to the original 
ession. 


he original expression is recorded 
itten). The digit representation of 
term is written, one term beneath 
other. The proper indications for 
dination or noncombination are in- 
din the following vertical representa- 
of the expression. Note that the 
al representation of the expression 
ms the same for both transforma- 
The only difference is that the 
ination and noncombination indica- 
are reversed. 


To AND Form 


ee DIC SBy A 
3 a ee 


Press 7 


-_~ 


(A) 


Fig. 1. 


(B) 


Sample expressions 


A—Relay diagram 
B—Static diagram 


To or Form 


Ia ID) (GIB Z| 
oy 2 Se 
+4444 
=f ea 


+—+—4-4+ 
gts ae eee 


+4444 
=e ea 


+—+—4-4-+ 


Rules for Transforming the Function 


1. Rules for any digit of two terms that 
are to be combined are: 


(a). AO in the digit of one term and a 1 
in the corresponding digit of the other term. 
The two terms cannot be combined and 
there will be no resulting term. 


Example: Digit 
1 
———=no resulting term 
0 
(6). A O in the digit of one term and a 0 


in the corresponding digit of the other term 
will yield a 0 in the same digit of the re- 
sulting term. 
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Example: Digit 


0 
ee) 


0 


(c). A O in the digit of one term and a 
dash (—) in the corresponding digit of the 
other term will yield a 0 in the same digit 
of the resulting term. 


Example: Digit 
0 
SEG) 
(d). A 1 in a digit of one term and a 1 in 


the corresponding digit of the other term 
will yield a 1 in the same digit of the 
resulting term. 


Example: Digit 


(e). A 1 in a digit of one term and a dash 
(—) in the corresponding digit of the other 
term will yield a 1 in the same digit of the 
resulting term. 


Example: Digit 


(f). A dash (—) in a digit of one term and 
a dash (—) in the corresponding digit of the 
other term will yield a dash (—) in the same 
digit of the resulting term. 


Example: Digit 
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2. Noncombinable terms are left in 
their existing form. 


3. Terms may be grouped that have 
consecutive noncombination indications 
between them. The terms will be re- 
peated (rewritten) and the indication 
between the terms will be omitted. It 
will be understood that the terms are not 
to be combined, 


4. The indicated combination concern- 
ing the smallest number of terms will be 
performed first. The successively larger 
combinations are then performed until 
no more remain. 


one or the other may be eliminated. 


Example: 


Gabe 
ned Wiens 
p= 


meh 


9. For two terms that have the same 
number of literals and the identical digit 
representation except for one digit. If 
for that digit one term has a 0 and the 
other a 1, the two terms may be expressed 
as one new term. This new term will 
duplicate the common digit representa- 
tion of the two orginal terms and a dash 


Sample: 
To AND Form 
POD. Ce BoA FAD“ Gap EsD CUB aA El DiGaPiaA 
eee pees oo ee etn Ha et ee rs hee 
ae a ae Wipe kee iy ae Toe er 
_ - 44-444 —— —-+—4—+—++—— 
tae st oe) eee Onan camel OAle= Ome 
——+— $$$ — 
OPP Ses he ieee: i ee Sexes) ain BigD GaskaA 
Liniget Steet 
See Pt Ney a tis ier 6 =. de 5 ae Ohh s Ole tah 
———_+—+—+_+_+—_— 
(ies a r= Lie ube has 
(Ee meh pe eal ve Mes eS 
Salyers ae. J 
ee ie 
+—+—+—4+—+ 
<b) VS Sa 
To or Form d 
eA COB BD Cr BeA EDUC EA EDCBA EDC BE 
2 eh Be te ler es or ee es ~ 
+—+—+—+—+ 2 
See ee en ae Sp eel ae ee | 3 Sa 
= Saas 2 SQUAT ELMS = ALR TT CR ———4+—+—+- +4 
eee = Ren ea jes pS, eee i Rae ete 0 1 2a 
+444 + 
~ Vi 2s ee Se So. vy SPREE oo Oe at ea 0 = 
+—+—+4+—4—+ 
a a --1-- --1-- 1-0-- 1-0 - 
ih eee Cie yes a = 34.0 ee 
++ +++ 
Sh ar Sue 1 ae hh a= i Near are Ap) 1 - =a 
+—+—+—+—+ 
ee = T ees UG eae ooh Dar pe ROD path WPA As - 1 = 
—+—+4+44— —+—+++—+4— 
22 ira ees att Wik wien aN ho wate ee een 255 ae 
Re ae ae 1 iT A eee 121 2 
ee Wa nee sey 
= hie a ry 


Rules for Simplifying the Functions 


A literal will be considered to be con- 
tained in a digit if that digit has a 0 or 1 
in it. The literal digit representation of 
a term is the 0 and/or 1 configuration 
contained in its digits. If a dash (-) is 
present in the digit the literal is not con- 
tained in that digit. 


1. If two identical terms are repeated, 
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(—) will be placed in the digit that varied 
from 0 to 1 (or 1 to 0). 


Example: 

EY DAG= BoA Eh DIG Brea 
en ne ee eo a | 
= 9 la So ot 


3. For two terms, one of which con- 
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tains a lesser number of literals than & 
other and also has its literal digit re 
sentation duplicated in the other te 
The term with the greater number ¢ 
literals is redundant and may be é x 
inated. 


Example: 
E 


E DG. Ba 
me ha Oa ee 

; 
4, For two terms, one of which cont i 
a lesser number of literals than , 
other and also has its literal digit re 


sentation duplicated except for ont 
in the other term, if for that ofie dig 
term has a 0 and the other a 1, the 
containing the least number of 1 
is retained. The other term bec 
new term that duplicates the co: 
digit representation of the two 
nal terms and a dash (-) is plac 


> digit that varied from a 0 to 1 (or 
0 0). 


ample: 

mac BA EDC B A 
foe LS a SL 
| [Uae =1-11 


(he sample for the anp form used in 
m 4 under “Rules for Transforming the 
action,” indicates that the remaining 
0 terms cannot be simplified and there- 
le constitute the final expression. 


fap CB A BUD 
= ae a toe 
See (4)0 1 
om 1 - - - (5) 0 = 
a (ay ie 
ma ~ 0 — - (7) 8 
mM 10 - - (8)(8)(9)(10) — - 
Se | i ge a 
mm 1 - — 1 G2)(18y 1 
ae 1 = 3 
or? 1+ = 
eet — = 
Se a 
EY Del Bua 
(2) a= sly = 
(4)0 1 ae 
(5)0 - - - 


i 
(61a Ou 
(3)(8)(9)(10) - - - 1 
eens) 1 = 


DUCTIONS IN THE BASIC PROCEDURES 
9F THE METHOD 


t is possible to shorten the basic 
sedure which has been outlined in 
This reduction is accomplished 
fetaining the basic rules of com- 
ation and simplification and by adding 
following modifications. 


In order to record the original 
ection for simplification and transfor- 


ion: 


The same procedure is followed as 
ined in items 1, 2, and 3 under ‘‘The 
ic Procedure of the Method” ex- 
that the digital representation of each 

instead of each product will be 
ten. 


iple: AB[CDE+ACE(C+D)] 
ap 0. BA BoD CBA 


en hie (io) 1 ee 
Me pDCBA EDCB A 
eee ot ie |) 


To record the function the digital 
‘esentation of each term is written, one 

beneath the other. A P will be placed 
t to all terms that are products. An 


ies) 
S 
Q 
yy 


Sample: A 
iT 
1 


|1orl 
ele l 
FOr | 
| 
nets 


(c). The logic functions between terms 
will be handled exactly as was outlined in 
items 4 through 7 under ‘‘The Basic Pro- 
cedure of the Method.” 


To or Form 


CBA E 
(2) — 
(4) 0 
(5) 0 
(6) 1 


bPely 


CO 
(3)(8)(9)(10) - 
ae (11)(12)(13) - 


Fay Ore ef 
| 
| 


Bere | OOF |] 
I 
aie 


E D 

a aah ie 

(5) 0 

(6) 1 - 

(3)(8)(9)(10) — 

(4)(7)(11) - 
(12)(13) 


| 
Wn aaicy 

! 

| 


— 
| 


Terms 3, 8, 9 and 10 are simplified by rule 3. 

Terms 12 and 13 are simplified by rule 1. 

Terms 11 and 12, and 18 are simplified by rule 3. 
Terms 7 and 11, 12, and 13 are simplified by rule 2. 
Terms 4 and 7, 11, 12, and 13 are simplified by rule 3. 


To AND Form 


EP D CBA 
es 


Sample: 


Tt 
Lee Le See 
a) SO VAL 
ae Mie eS 
To or Form 


EDGES B eA 
COA Tz 


ee ee ee ee oe 
IF Melee 5eP 


= Ll a! eS 


2. Combination. The basic rules are 
applied for indicated combinations with 
these additions for transformation to the 
AND form: (a). A term that is a product 
will be combined as-a whole term; (6). 
A term that is a sum (S) will be combined 
digit by digit. 

For transformation to the or form: 
(a). A term that is a sum will be com- 
bined as a whole term; (6). A term 
that is a product (P) will be combined 
digit by digit. 
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To illustrate the combination consider 
two products 


1 AID ATORS Toh. 92) 
DBAS—) Vea ee 
ED 0b = shee 
and two sums 

IP IDG 18h VAN 
(E+B+4)=0 - - 10 S 
(D+ CHB) =o 0 one 


Transformation to the sum of products 
(AND) form: 


Example 1. Combination of two products 
ESD CBA ED Germ 
Se Ay oul CAP) Oeil ee ae 
0 1- 1- P 


Example 2. Combination of two suris 


Jie JOY ENED Al BD Ca Bupa: 
0 = = 140. San 0 Oe 
oa ene OF 
= 01 SOS emp eares 
22 ()) Sean ees 
ese teas SP ee 
Bre rey (get | 
40h Sane 
2 5 gil aan) 
oy etn) 
Example 3. Combination of a product 
AND sum 
Belo CoB aA EID NG Ae 
a Sar Stile G2 COPS bee al 
OS LORS Se 
Transformation to the product of sum 


(or) form: 


Example 4. Combination of two products 
EOD CBA iD Ce Darl 


aad) Se eel ee 0o1l1-- - 

OMA = Ae ee - l1-- = 
0--i1- 
Nn eee 
Ee ee 
0--- 1 
- 1-- 1 
- - - 1 ji 

Example 5. Combination of two sums 

‘By ID Gab eA BOS ON (OM ey 71 

Oo Fate Ae Oi AS: OM On ie tes) 


ONL aera 


Example 6. Combination of a preduct or 
a sum OR a product 


BoD Ge BIA 


ED ERG: 
eo idle en Mallya ot Oa 


_ 


B 
Ou 
0 


Oye OP = al 


Fol Gas tg 
=) 


O12 = 


— 
| 
a) 


1 
BA 
0 


SISES City So 
Bhyr 
Ee 
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Upon close inspection of these examples it appears that some of the resulting terms 
‘of an indicated combination do not appear. This is due to part (a), item / under 
“Rules for Transforming the Function.” 


3. Simplification. The examples in the preceding section were not simplified. 
Further inspection of example 2 yields that one term contains only the B signal. 


By rule 3, all other terms that contain B (a 1 in the B digit) are redundant and 


may be eliminated. This simplification may be accomplished at the same time 
the combination is made, thereby directly reducing the number of resulting terms. 
To accomplish this, the following rule may be used: 


If in performing the indicated combination (for either transformation) of literals or 
whole terms, one of the items has its digital representation duplicated in the resulting 
term, all other combinations of this item will prove redundant and therefore need 
not be considered. 


Example (from Example 2) 


PSS) 
fee le 
! 
i 


Example (from Example 6) 
2, AD (Gey 2 JE IBY ME: Ves IA DG BrAl 
OY Cai 0) Ih al @ 


(yee a PS 1 PS 68) - (2)(3)(4)0 . 
(3)0 — 1.0 5 OS" O41 0-101 
(4)0 1 --1-P (2 re on cat UN nei? 
Sample: 
To AND Form 
(1) (2) (3) 
BEE GiB eA BED SCN Bea B-D-C-B A 
(2)--- - 1 1P a ae oem | UTS heehee 
(CD i ne a ete ay (B) dee SSA ep. Oi teO er 
—+—4+—4-44— 
(COO ple PO (Ga) ah On ee al 
(Olt Vhs elas S 
To or Form 
(1) (2) (3) 
ED) Cee BiaeA BDAC Abad EDS CoB A 
(Fea eee Bee bn (D\e i SC ep a@a----1 
———+ —+—4-44-+ b= = = 1.5 
(Byrael at - Pp 1 = 0% = CS ORE 
Es = 
(4/0. =-.0-=-1 P 1 See =o SANS VOC ewe =) 
+—+—+—-4-+ f=10-- 
(25) CS ed te AE ai eo g-1--1 
(3)(4)(5)- 1 0 - - b0Q=-t2 2 
=" LP MS (i ae 
0 - Le=- = Geoy: (i ae 
sri ali ee: 
=p et nl se 
(4) (5) (6) 
Ee DSC BivA DUCE aA ED CBA 
be GL @a- -- = |] C=! 
b- - - 1 - ee ha  tpbes cel yet aed 
cl - O - = el = 10 = 25 ol = 0) = = 
302169105) AUR a a ae ae RO! (eae 
(et 0 ea fi Os ea eee Si a ean 
bhO - 1 - - fg=-1 = = 
g- 11- - 
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OBTAINING THE MrinrmuM Form OF A 
SOLUTION 


After all indicated combinations 
been performed the resulting t 
constitute a solution. This solution 7 
not necessarily in the minimum 
and the following process may be u 
to obtain a minimum form of a solutior 
The process will indicate that a term 
redundant if its function is provided 
in two other terms. In order for t 
to be valid the two terms must cont 
fewer or an equal amount of literals th 
the third term. The process is used 
follows: 


1. Consider the terms containing 
greatest number of literals first. 


(a). Substitute a 0 and a 1 in a ¢ 
containing a dash (—) in the considered te 
This will yield two new forms of the « 
sidered term. 


(b). Compare one of these new for 
the remaining terms. If it duplicates 
literal digit representation of another t 
the other new form of the considered t 
will be compared. If this form dupli 
the literal digit representation of anot 
term the considered term is redundaa 
will be eliminated. 


(c). Should the first 0,1 substitution # 
prove the term redundant the proceé 
will be repeated for another digit con 
a dash in the considered term. Th 

continue until the term is determi 
redundant or that it is to be maintaine 


2. The next term containing 
greatest number of literals is then 
sidered and checked for redundancy 
the same manner. This is done t 
there are no more terms to give 
necessary comparison, i.e., that there 
not three terms that have at least 
equal amount of literals. 

The sample function for both tr 
formations is already in the minir 
form. 

The following examples will iltust 
the procedure. 


Example 7: 2 
Cl Bea C Be 
(2) - 1 1 (2) -> i 
(3) ee= ao (3) 1 S308 
S44) ae " 


All terms have the same num be 
literals so that each will be considere 


Term (2) 
Making the 0,1 substitution in 
digit yields the new forms of 0 11 
desi 


Comparing 0 1 1 to terms (3) an 
shows that neither term has its 1 
digit representation duplicated, 

fore term (2) is not redundant and 
be retained, 


m (3) 
; aking the 0,1 substitution in the B 


omparing the 1 1 0 form with terms 
) and (4) it duplicates the literal 
igit representation of term (4). The 
| 1 0 form is then compared to term 
2) and it does not duplicate this 
tm’s literal digit representation. 
herefore term (3) is not redundant 
d must be retained. 


n (4) 
aking the 0,1 substitution in the A 
git yields the new forms 1 1 0 and 
iv oe 


omparing the 1 1 0 form to terms (2) 
d (3) it duplicates the literal digit 
presentation of term (3). The 11 1 
rm is then compared to term (2) and 
duplicates the literal digit representa- 
on of term (2). Therefore term (4) 
redundant and is eliminated. 


ple 8: 

DC B A & D BA 
Oo =O 1} (2) = 40) = “01 
ete i (ay 2 4 a 1K 
Peete (5) =O ai = S&S 
eee) a(t) 0" 1 Oo) = = 
fa- 0 Q 1 

010- - 


ibstitution of a 0, 1 in the D digit 
erm (6) yields the new forms which 
icate the literal digit representation 

(2) and (7). Therefore term 
$s redundant and is eliminated. 


(A) 


(B) 


Substitution of 0,1 in the D digit of 
term (3) yields new forms which dupli- 
cate the literal digit representation of 
terms (4) and (5). Therefore term (3) 
is redundant and is eliminated. 

The other terms are checked for re- 
dundancies and none are found. There- 
fore all must be retained. 

After becoming familiar with this 
process the substitution and comparisons 
will be made mentally and the redundant 
terms eliminated by inspection. 

This process allows a term that is pro- 
vided in two other terms to be eliminated. 
Elaboration of the process so that a 
function that is provided in three or 
more terms can be eliminated. Generally 
the elimination of these redundant terms 
will yield the minimum form of the 
solution. Many other methods (includ- 
ing those outlined in the Karnaugh, 
Quine, and McCluskey publications) may 
be used to obtain the minimum form of 
the solution. 


WRITING THE FINAL EXPRESSION 


After all indicated combinations have 
been performed and all simplifications 
have been made the remaining terms are 
in a group. For the conversion to the 
AND form each remaining term will be 
a product and the terms will be separated 
by or functions. For the conversion 
to the or form each remaining term will 
be a sum and the terms will be separated 
by AND functions. 

For either conversion each term will 
be written as follows: 


moow,r 


eo 


mMmoowr 


Ol 


mo 
. (3) omD> 


Fig. 2. Relay and static diagrams of transformation of sample expression 


A—AND form 
B—OR form 
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1. The variable will be present if a 1 
appears in its digit of the particular term. 


2. The negation of the variable will be 
written if a 0 appears in its digit of the 
particular term. 


3. The variable will not be written for a 
particular term if a dash (—) appears in its 
digit. 


Sample 

To AND Form 
ie) Cae Bae, 
ls PR A ee Arche CD 
OF a OMe rr ABE Gs Di 


Final expression: X =ABCDE+ABCDE 


To or Form 


CB eA 
sacl B 

i (C+E) 
0 


(C+E) 
ort A 


i Geet is 


X =B(C+E)(C-+E)AD 
=ABD(C+£)(C+E) 


Final expression: 


Figs. 2(A) and (B) show relay and 
static logic diagrams of the sample ex- 
pression for both transformations. 


Summary and Conclusions 


In summation, two new methods of 
simplification and transformation of Boo- 
lean functions have been presented in 
this paper. These methods are utilized 
to obtain a minimum solution for a 
desired form and may be referred to as 
conversions of the original function. One 
method is used to convert the original 
function to a minimum sum of products 
and the other method is used to convert 
the original function to a minimum 
product of sums. The two methods 
parallel each other in that they utilize 
the same basic rules, the major difference 
being the recording of the original 
expression. Because of this similarity 
the two methods have been explained — 
as one basic procedure. The funda- 
mental principles of this basic procedure 
have been outlined in detail in the discus- 
sion. A reduction has also been ex- 
plained that allows the basic procedure 
to be considerably shortened. A 
“sample” Boolean function and inde- 
pendent “‘examples’’ have been included 
to better explain the procedure. Examples 
of conversion to both forms of two 
Boolean functions have also been in- 
cluded in the discussion. 

The following advantages occur 
through the use of the procedure pre- 
sented in this paper: 


1. The method for either conversion is 
based upon the postulates and theorems of 
Boolean algebra. Therefore, (a) any exten- 
sion of the procedure or any portion of a 
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similar process that is also based on Boolean 
algebra may be incorporated to expedite 
the procedure; and (b) the manipulation of 
the original expression in lettered variable 
form by standard Boolean algebra methods 
is eliminated. 


2. The original Boolean function appearing 
in either pure form (sum of products or 
product of sums) or a combination of both 
forms may be converted directly. There- 
fore, (a) it is not necessary that the original 
function be initially presented in either 
pure form to accomplish minimization or 
conversion, as it is in the Karnaugh map 
method; and (0) it is not necessary to 
analyze the original expression for the basic 
transmission functions as it is in Mc- 
Cluskey’s method. This eliminates 
duplication of information stated initially 
in the original function. This duplication 
is encountered in the combining of these 
basic transmission functions to achieve the 
minimum form of a solution. For either 
conversion by the new methods the original 
function is the true base. 


8. No restriction is placed on the number 
of variables or the complexity of the 
Boolean function for either conversion. 


4. The basic procedure is the same for 
either conversion. Therefore by becoming 
acquainted with one of the new methods 
of conversion it can be applied easily to 
accomplish the other, eliminating the 
necessity of learning two distinctly different 
methods. 


5. Although not a new concept the 1,0 
and dash (—) variable representation 
utilized presents first, a clear, concise and 
easily handled representation of the original 
function for either conversion; second, a 
step-by-step tabular form for either con- 

- version that will reduce the possibility of 
omitting literals of terms or whole terms; 
and third, a record that can be checked 
conveniently for errors. 


6. The procedure is such that for either 
conversion most of the redundant terms 
are disregarded or eliminated as the con- 
version progresses. This means a minimum 
number of redundant terms are carried 
through to the last step. 


In conclusion, from the fundamental 
precepts of the procedure presented: 


1. New parallel methods of simplifying 
(obtaining a minimum form) Boolean 
functions have been provided for static 
and conventional relay-circuitry reduction. 
The equipment itself determines the con- 
version necessary to provide for maximum 
savings. 


2. The use of standard Boolean algebra 
operation on the function has been elimi- 
nated in these new methods. Some of the 
restrictions encountered in using other 
methods (previously mentioned) have also 
been eliminated. 


3. The rigidly systematic parallel methods 
and tabular 1,0 and dash (-) variable 
representation present an easily pro- 
grammed process for digital computer 
application. 


Extension of these basic ideas will 
surely provide even further reduction of 
the procedure. 
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Appendt: 


The preceding sample and examples considered no more than five variables. 
following examples illustrate how both transformations may be accomplished. for express 
containing more than five variables. The method utilized in these examples is the reduct 


of the basic procedure. 


Example 1: (C+D+BF)[AB(C+E)+BD]+ABC((D+E+F)+BF] 


To AND Form 


FCRE DEB tA F. BI DVC~ BA P ECD COBSA 
ver lt BAER ---1-- ---1- 5. 
+—4+—4+—+4—4— + 
~ = 0f- — - Pp --0-- - --0-- - 
4+—+—+—+—+—+ 
POS PN 1 = 4-1 - 1---1- 
yah ses oL ai Ba uae tiene ee 
St Pee ee eS == = Oe tet - +-- 0141 
+4+—+— +--+ 4+ 
ee. Ey het. COME AR --1-0- --1-0- 
$444 | 
ees Ue cO al te 1 =h Ot ---+10 
PONG hic ongl got es =" 0-1 0es 29 Peat Cee 
+—4+—4+—44—+ OO ae ae 
0---0- 2 be Tate ad = Uh «ld /iee 
PORNO Bae A PLES DCs FE DC B Ave 
eet te se rie 1 al =. 2 So ee =- 1.- 1,251 
cot ne peli janes shale of uaa -~ = 4 1). 
aed Verio ele ~ = O01: 28 == 0 O L 1g 
Sea Opal and 1 = SO at 1 ~ 0.1 23 
tae era erie el - i hia = ee 1 0 1 
1 TOAST t 
Ste aero yO 
X=ACE+BCD+ABCD+ABCF+ABC 
To or Form 
ToRW ay spre tg Ry RoE D SG) Boe FE D CBA 
SNeet= shou | eas 1 = 0. PS 1-o1- - 
= S00) s-¥ Seg == 01 We -— =. Pas 
Pee ace) teal, P = = “hee es ee 
—+—+—4+— 44 4—— 
----1i1 P =) Sid's a ee a ee | 
+—4+—4++4—4—+ 
SO se Othe! Sam ES = Dial 2. She Oe 
Saha ey =P ae bbe gy) Sion = = 11 0p ee 
=<"“1°= 0 0 = = 1 0:0 = 
fae te = ee Re P 
—4+—+-4-4-4-4-— ---1-- -=--+1+8 
1 O74 Ate BAG = 2a ee, eg 
- = ‘— 1 ees - 1 
- ———— 1 “Ose” = Fone 
OG vee Nig) ae, ae 
ey ROC VE PRD" CAR A 2, 
Ts S0p alipe= a 1 - 0 1 - - (C+D+F) 
Hy ee --- 11 - (B+6) 
----+= 01 - - + 0/1 G+ 
= ees ~ L040 =X Beis 
a <a += Mosabe angae es 
ee ear Bie a te ( ) 
Sib LO 0, tee 
ioe Ol.” ee ee 
and. 1 Sor wa 


X=(C+D+F)(B+C)(A+B)(B+C+E)(A+D) 


Example 2: (The digits representing the variables will be reversed to those in the pre 
cases. This has been done deliberately to illustrate that the variables may be ini 
assigned to digits in any order. By assigning the corresponding variables to the ec 
digits throughout the procedure the order has no effect on the conversion.): 3 


= 


X=[ABC(D+E+F)\(C+G+H)|DG+[AGH+(A4+B+F\(C+D+GH)] — 
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Discussion 


Peter L. Randlev (IBM Research, York- 
town Heights, N. Y.): While Mr. Howard 
has been correct in attempting to find a 
simple straightforward method to obtain 
the most economical switching circuit, his 
method does not consistently do this. 

The most economical form of a combina- 
tional circuit is dependent upon the com- 
ponents to be utilized. However, in de- 
veloping a method for determining the 
configuration of such a circuit, it is not 
sufficient just to specify whether an alter- 
national (sum of products) or a conjunctive 
(product of sums) form of logic shall be 
used. Often, the most economical con- 
figuration should be a combination of both 
these forms. It is also important to make 
provision for additional rules to obtain 
this configuration in view of: 1. the 
cost of each term in an expression as well 
as the cost of each literal, 2. the rela- 
tive cost of terms composed of various 
quantities of literals when, as is most often 
the case, this cost is not linear with the 
quantity of literals, and 3. the cost of a 
circuit configuration composed of several 
levels of logic as would be obtained if the 
function were in a factored rather than a 
normal (pure) form. 

Any method that is to be able to find the 
most economical circuit configuration must 
be able to take into account factors such as 
these for defining a minimum. The author 
has no definition of a minimum. 

For any method to be easily and success- 
fully programmed for a digital computer 
application, as Mr. Howard suggests, it 
must be both systematic and unambiguous. 
As presented, his method fails this criterion. 
The rules are not stated explicitly enough 
and, rather than being rigidly systematic, 
the examples themselves are demonstrations 
that the rules allow any of several sequences 
in performing the operations. To be used 
with a digital computer, it would obviously 
be necessary to use some means other than 
lines of several lengths to indicate the 
sequence of groupings required in trans- 
formation. This system does not even 
appear to be the best with hand methods 
as the comparison of the relative lengths of 
the connecting lines would become quite 
difficult to see in large problems. However, 
the author should be able to modify quite 
simply his rules and notation for use with a 
digital computer. 

Other than the purely clerical operations 
of recording the original function and 
writing the final expression, the method 
may be divided into three general opera- 
tions: 1. transformation or combination, 
2. simplification, and 3. minimization. 
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The transformation operation changes 
the original function into a normal (pure) 
form. As desired, this would be either 
the alternational normal form or the 
conjunctive normal form. The simplifica- 
tion operation apparently is intended to 
find the terms that have the fewest possible 
literals. This is accomplished by pairing 
the given terms and removing those terms 
that are found to be completely contained 
within one other term. However, the rules 
do not always allow the term with the 
fewest literals to be found. The following 
is such an example: 


Example A: 
Original data Desired answer 
Ome 0) - 01 0 
=o lO 21 =F! (Oo ot 
Sle ol 1-- - 
1 - 0 0 
ig Pale ile 
0 0s — 


The stated rules will not allow any simpli- 
fication of this original function. 

It is probably Mr. Howard’s intention 
to define minimization as that expression 
of the function whose sum of the total 
quantity of literals in all terms plus the 
quantity of terms is the smallest possible. 
As the author has apparently recognized 
(in the last paragraph under “Obtaining 
the Minimum Form of a Solution’’), his 
method does not necessarily find a minimum 
solution under this definition. The opera- 
tion the author has described for use in 
minimization is to remove redundant terms; 
that is, terms that are expressed by two 
other terms. This method often removes 
a term that should appear in the final 
answer. For example, see the following: 


Example B: 


Original data Desired answer 


i SO) i = ih OP Wee = 
re Ogee el leesl@e Al 
IL al a 30) 1 Onn: 
OPS Ate al (Oy sedi ak 0) 
OO) = ah a OO; = i Bl 
I SO 

= (0) 2 ah = 


If the rules are followed consistently, the 
first term of the original data will be re- 
moved by the last two terms. which will 
themselves remain intact. This non- 
minimal answer would be 


1 


lrReoorl 
1Lolrr 
feoocr 


el ee ee | 
KFPOrr OO 


0 


In other cases these same rules fail to 
remove redundant terms, as in the following 
Example C: 


Original data Desired answer 


jeer ihe ios 
powered 1 = OO 
oN oSleiat = 1 Sieeo 
a oP eek = nr. ol 
{one 


When the rules are more precisely stated 
as to the sequence of operations for use with 
a digital computer, the sequence of the 


original data becomes important in 
tempting to find a minimum soluti 
This is demonstrated by example D. 
this problem is “minimized” with the d 
both in its original and in an inverte 
sequence, two answers of different c 
are obtained. } 


Example D: 7 
Nonminim al | 
Original Desired “minimized’ 
Data answer answer 
tee 1 1 = 1 1, 20.32 
t-0 hb =A1 0 We =]. 0s 
Ol OG -=4 0 +10 == alee 
Ini] sie - 10 
- 101 


In the reference 2 cited by the autho 
Professor Quine discusses the minimizat 
of problems he calls “‘irredundant.” Pro- 
fessor Quine’s example is ‘ 


Example E: 


Original data Desired answers 


i = 1.0 = S6.0me 
jee o - i or iam 
Ee 0) ~*1°°0 ~ Oe 
heel) We 


As Mr. Howard’s method is comple 
dependent on finding redundant terms, 
cannot find the minimum solution fr 
these irredundant functions. The desi 
answers of these five examples can all 
checked by other means (such as a tre 
table) to insure that they are equival 
to the original data. These problems 
only samples of those problems that cai 
not be minimized by the author’s rules. ~ 
It is my experience that as the problem 
become larger and more complex, finding: 
a problem belonging to a class which canno 
be minimized by this method becomes a3 
more frequent occurrence. For a fr 
minimization process to be useful, it m 
be possible to recognize the nonmin 
cases even if a minimum answer cannot al. 
ways be found (so that they may be solved 
by other methods). 
There is no step in the author’s method 
in which the intermediate results beco 
more complex in form than the origina 
Thus, rather than referring to this as 
minimization process, it would be more 
acceptable to call this portion of the author's 
paper a method to remove redundancy, 
Extensions of the method to allow 
first, to recognize cases where it cannot fin 
a minimum, or second, to find a minimu 
answer for all cases, will make the proce 
more complex and lengthy rather the 
reduce the necessary procedure. It 
also necessary that a minimization pr 
cedure include the ability to handle “d 
care” conditions. These are norm 
created by specific input combinations th 
can’t exist and taking them into consi 
tion often allows drastic reductions in 
required circuitry. 
This method becomes, in common wil 
all other methods known to this discu 
extremely difficult to use by hand f 
problems much larger than those thi 
can be simplified by inspection, just bec 
the problems are large. However, 
large-problem area is the one in which f 
greatest savings may be made by the wi 
of simplification schemes and whic 
therefore, requires the most work. 


Jury 19% 


lw. Stuart (Industrial Nucleonics Cor- 
ration, Columbus, Ohio): Mr. Howard 
6 done a thorough job in presenting a 
thod of converting and simplifying 
binational logic circuits. Circuit de- 
a based upon the use of logic functions 
B received increasing attention in recent 
mrs. 

1 he binary nature of the normal industrial 
ptrol provides a comparatively simple 
tical system. It can be shown that the 
s of ‘symbolic logic” are applicable in a 
her restricted sense. Specifically, de- 
es capable of providing the logical 
bperties of assertion and negation are 
p basic requirement. Assertion corre- 
bunds with on while negation is representa- 
e of OFF, 

Assertion can be categorized into two 
cipal groups. Alternative assertion is 
scribed by the verbalization or whereas 
ordinate assertion (coincidence) is de- 
iptive of AND. Negation is by definition 
Overt representation of the absence of 
ertion. This inversion is often termed 
e complementary function. The prac- 
al importance of grasping the concept 
negation cannot be overemphasized. 
gation is not a passive zero state. For 
ample, a normally closed contact does 
provide negation unless it is also con- 
ted to a voltage source. Furthermore 
unoperated normally open contact is 
synonymous with negation, only com- 
tible with it. An active signal must be 
bvided to represent the open condition 
this normally open contact in its un- 
stated state. Identical reasoning applies 
static switching devices. 

t is possible and practical to construct 
ire logical systems, including memory, 
lizing only negation and one of the two 
ettions. Therefore, a given system 
ght contain only or and nor. An 
ertnative system could contain only 
ib and not. It should be evident that 
system lacking negation, or NOT, is not 
sible. 

system composed of all three basic 
ie functions is defined as an English 
ic type. It provides the most flexible 
stem design choice and is apparently the 
sest approximation of normal thought 
evidenced by the language structure of 
least the English-speaking people. From 
practical standpoint, such a logical system 
easier to understand, design, and trouble- 
ot. This consideration is important 
ce the basic description of the operation 
quence is invariably presented in a 
alized form. Furthermore, the usual 
atrol design methods depend on this 
ibility unless a mathematical tool such 
the algebra of logic (Boolean algebra) 
itilized. Thorough mastery of this tool 
necessary to allow effective design. 
is algebra is limited in over-all applica- 
n since time variant or sequential opera- 
n is difficult to take into account. There- 
e, the entire class of memory types is 
uded. However, for optimizing ex- 
mely large systems or very complex 
ings of combinational logic configura- 
ns, this algebra is very effective and 
rth while. 

The usual starting place in tackling a 
itrol design is with a formulation of the 
iditions necessary to activate each 
put. Upon completion of the individual 
cuits as well as after final integration of 
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the over-all system, it is desirable to expend 
a reasonable effort in checking and simpli- 
fication. The art of simplification is tied 
closely to the economics of alternative 
logic components available. In the case 
of the static-switching system under con- 
sideration, the most economical system is 
generally achieved by the maximum use 
of or units. Furthermore, the manner in 
which memory units are formed has an 
important bearing upon formation of 
optimum circuit design. Broadly, two 
processes can be used to affect simplification 
of combination circuits composed of ors, 
ANDs, and NOTs or variations thereof. 
These processes are conversion and in- 
version. Conversion is concerned prin- 
cipally with OR-AND circuits, while inversion 
is concerned with nor functions as the 
term implies. 

A given portion of a combinational 
circuit may be the termination of two or 
more ANDs at a single or. Conversely, 
it may be the termination of two or more 
ORS at a single AND. It can be shown that 
a one-to-one correspondence exists between 
these alternative circuits such that one may 
be converted to the other with about an 
even chance that simplification will result. 

The greatest possibility for worth-while 
simplification exists whenever a _ given 
combinational circuit contains a repetition 
of one or more identical signals. The 
final circuit may or may not contain re- 
petitive signals since economy is not 
primarily a logical consideration. It is 
equally important to note that if a signal 
A and its negation or complement A when 
found in a given combination are often but 
not invariably redundant. 

Frequently, the availability of the 
normally closed contacts of sensing devices 
permits simplification by permitting inver- 
sion of the circuit. Inversion entails 
changing OR conditions to ANDs and AND 
conditions to ors with all input signals 
inverted (negated) and providing the 
complement, or NOT, of the resultant circuit. 

Both conversion and inversion may be 
carried out on a given circuit for the 
ultimate in simplicity. The type of memory 
circuit with which the combinational 
circuit is connected in a large measure 
determines whether conversion, inversion 
or both are profitable. 

This paper presents a complete and 


usable method of handling OR-AND con- ~ 


versions and performing a straightforward 
procedure of simplification. 

Memory circuits would appear to be 
most effectively handled with sequence or 
time-bar charts to define the combinational 
functions which occur momentarily that 
must be remembered for subsequent use in 
the circuit operation. Unless such a 
method is used to reduce the theoretical 
memory possibilities to actual require- 
ments, the task is almost unmanageable. 
At any rate I have not seen a really practical 
method of handling memory functions with 
the use of Boolean algebra alone. The 
analysis of memory situations certainly 
should present a fruitful field for future 
investigation. 


R. L. Howard: The author wishes to 
express his gratitude to Mr. Stuart and 
Mr. Randlev for their discussions. 

The author’s background is basically 
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that of industrial-control systems while 
Mr. Randley’s is predominantly concerned 
with digital computer applications. His 
discussion therefore provides an appraisal 
of the procedure from the computer view- 
point, which is an area with which the 
author is least acquainted. Primarily Mr. 
Randlev’s comments are directed to certain 
generalities contained in the paper that 
would seem to imply that the scope of the 
procedure presented is unlimited. It was 
not the author’s intention to create this 
false impression, particularly regarding the 
minimization of Boolean functions, 

The procedure presented does not provide 
methods for obtaining the most economic 
solution for a particular application or 
function. Surely many factors must be 
taken into consideration before the optimum 
form can be determined and it has been 
assumed that the designer has solved this 
portion of the problem before utilizing the 
procedure. Those applications that do 
find their most economic solution in one 
of the pure logic forms can be resolved 
directly from the procedure. Other applica- 
tions that find their optimum solution in 
the combinational form, will require addi- 
tional effort on the part of the designer 
in using the procedure as is the case with 
all similar methods known to the author. 

Different lengths of lines are utilized 
throughout the procedure to indicate 
groupings within the function as well as 
serving as operational devices. This does 
not mean that the procedure is restricted 
to these symbols any more than the AND 
or or functions in Boolean algebra are 
always represented by the same symbol. 
The author readily agrees that other desig- 
nations can be utilized for the desired 
indication and that these may allow easier 
handling of the procedure. 

The ‘“‘don’t-care’’ type of term has not 
been included in the procedure and what 
would be involved in taking this into 
consideration can only be determined by 
further investigation. In the same vein, 
it is the author’s contention that only 
further examination can determine the 
length of the procedure necessary to enable 
it to derive an absolute minimum for a 
particular function. 

It is a well-known fact that vast amounts 
of study have been directed toward the 
area of minimizing Boolean functions. 
Only a fraction of this effort could be 
covered in a paper of this length. Being 
aware of this, I made an attempt to qualify 
the rules for minimization, but obviously 
this qualification was not stressed to a 
great enough degree. Mr. Randlev’s com- 
ments on the simplification and minimiza- 
tion portion of the paper have pinpointed 
the restriction of the rules covering these 
sections and have brought these limitations 
firmly to the attention of the author and the 
reader. These comments (including the 
author’s presumed definition of minimiza- 
tion) are accurate and are invaluable to 
the author and the paper. 

In summarizing Mr. Randley’s discussion, 
I do not feel that the rules presented in the 
paper have been shown to be incorrect but 
rather that they do not include all cases. 

Mr. Randlev’s and Mr. Stuart’s discus- 
sions have been of great assistance. It is 
constructive analysis such as this that 
will increase the value of the procedure 
as a tool in logic design. 
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A Stability Grenen for Nonlinear 
Systems . 


Y. H. KU 


FELLOW AIEE 


N ANOTHER PAPER! Wolf reports 

on the recurrence relations in the 
solution of a certain class of nonlinear 
systems. In Wolf’s dissertation? a mathe- 
matical theory for the analysis of a 
class of nonlinear systems has been de- 
veloped. This paper reports on a sta- 
bility criterion based on the afore- 
mentioned dissertation and its application 
to nonlinear as well as linear feedback 
control systems. 

The method of stability analysis con- 
sists essentially of first finding the re- 
sponse of a nonlinear system to a given 
forcing function by either the method of 
partition given in reference 1 or the 
Taylor-Cauchy transform method.* 
These methods give the solution of cer- 
tain class of nonlinear systems!? in' the 
form of a power series, such that the 
response «(#) is governed by a function of 
the exponential type given in the follow- 
‘ing equation. 


x(t)= > hint” (1) 
n=0 


where h,, are constants for different values 
of n. 
The Laplace transform of x(t) is 


x(s)=)) (2) 


where s is a complex variable. 


Stability Criterion 


The stability criterion is given as 
follows. Given «(#) the solution of a non- 
linear differential equation (which repre- 
sents a given nonlinear physical system) 
such that x(#)=0 (e”). Then x(#)—0 as 
to if and only if the convex sin- 
gularity hull (CSH) of X(s) is entirely in 
the left-half plane (LHP). 

The CSH is constructed by joining the 
points of singularity of the transform 
X(s) inthe s-plane. It may be noted that 
the singularities of X(s) are not necessarily 
poles. 

The nonlinear differential equation can 
be given in the following form 


Z(D)x(t)+ F(x, x’, ...)=r(t) (3) 


where 
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A. A. WOLF 


MEMBER AIEE 


Z(D)=linear differential operator 
F(x,x',...)=nonlinear function of x, its 
derivative x’, ... ’ 

r(#) =the forcing function or input 

x(t) =c(#)=the response 

D=d/dt=the differential operator 

D?=d?/di?, ... 

x’ =dx/dt=the derivative of the response 
with respect to time ¢ 


A nonlinear system is said to be stable 
if to every bounded-decaying driving func- 
tion or input and for all initial conditions 
the response x(#) approaches zero as time 
increases to infinity. If the foregoing 
condition is not satisfied, the system may 
be oscillating in the steady-state or may 
run away and become unstable. 

In a linear system no mention is usually 
made of bounded values when speaking of 
the stability of the system. The reason 
for this is that a linear system has a 
response whose form is essentially in- 
dependent of the initial conditions. In 
a nonlinear systein, it may go from a stable 
to an unstable state during some transi- 
tion, such as switching which may alter 
the initial conditions. In such cases the 
nonlinear system is considered condi- 
tionally stable and the range of input func- 
tions and initial conditions must be 
specified for stable operation. 

In linear systems the CSH is composed 
of poles for lumped parameter systems. 
For stability these poles must be con- 
fined to the LHP. In nonlinear systems 
the CSH is composed of singularities 
(which are not necessarily poles) and for 
stability these singularities must be con- 
fined to the LHP. Thus the stability 
criterion given previously holds for non- 
linear as well as linear systems. 


Construction of the CSH 


The CSH and other related concepts 
in this paper are obtained from entire 
function theory.4 For the purpose of this 
paper, however, the following simple ex- 
amples show how CSH can be obtained for 
stability studies of physical systems. 


Example 1 
—. 1 

X(s)=)) ya (4) 
n=0 


This is a special form of equation 2, with 
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_ where a is a constant. 


han'!=1. Equation A or 2 (in general 
looks like a Laurent series. Ho 
one must not conclude that at s=0, . 
has an essential singularity. Instea 
may be that X (9 at s=0 is a regula 
point. 

Clearly equation 4 is an expansion of 


1 

X(8)= = 

At s=0, X(s) =—1, while at s=1, X(¢ 
has a pole. So the series of equa 
converges outside of the circle of con 
ence whose radius is R=1. In ge 
if the power series of equation 2 con 
outside of the circle of convergence whos 
radius is R=a, one must not conel | 
that at s=0, X(s) has an essential singiz 
larity. 

Since the singularity s=1 is not in LE 
the system is unstable. This can 
verified by taking the inverse transfor 
of equation 5. 


Example 2 


fee) a” 
X(s)=)) aa 
n=0 


Compared with equation 2, hyn! 
Equation 6 is a 
expansion of 


\ 1 
X(s)=— x 
orihe 
At s=a, X(s) has a singularity, which is§ 
pole. The system is unstable. 


Example 3 

il ua n 
X(6)= >») le 
is an expansion of 


af 
So Ra 


At s=—a, X(s) has a pole, which 
LHP. Hence the system is stable. 


Example 4 


. Let the simplified from of a given po 
series be given by ; 
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Im s= jw 


1 1 
ie s42' 31-3 


1 
Secepqaews) 0) 
e singularities are at s=—2, s=—1+ 
and s=—1—j2. Connecting all three 
mts gives the CSH in the form of a tri- 
sle. Since the CSH lies entirely in the 
P, the system is stable. 

he foregoing simple examples are 
en from the linear systems so that the 
struction of the CSH can be demon- 
ated in a straightforward manner. 
. 1 shows a CSH for a stable system, 
le Fig. 2 shows a CSH for an unstable 
tem. 


dius of Convergence and 
Direction of Singularity 


‘rom equations 4, 6, and 8 in examples 
}, and 3 of the foregoing, it is seen that 
Origin in the s-plane is not included as 
ngularity. In example 1, the circle of 
vergence has a radius R=1. In ex- 
ples 2 and 3, the circle of convergence 
a radius R=a. It is noted that in 
mples 1, 2, and 3 the singularity lies on 
circle of convergence. In general, the 
ius of the circle of convergence is given 


(11) 


lim |! hn|/” 
M—>o 


xample 1,equation 11 givesR=1. In 
mples 2 and 3 equation 11 gives 


fim | 11 hn|/"= tim |a"|/"=0 (12) 
I—> co m—>o 


direction of the singularity (which lies 
he circle of convergence) is denoted by 
angle 6, which can be obtained as 
ws. Let 


Rk 
5 (2) 


%=0 


(13) 
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A 


rR S- plane 


Fig. 1 (left). 
CSH for a stable 
system 
x 
ee 
orm 
Re seq- 
Fig. 2 (right). 
CSH for an un- 
stable system. 
Set of points K 
and its projection 
on 0=arg s 
where 
n! hy 
Bati= Pati 
and 


R\ Ok! nb 

n} ni(k—n)! 
Then 6 is the direction of the singularity 
with radius R if 
jim |x| ">2 (15) 
With the radius of the circle of conver- 
gence obtained from equation 11 and 
with the direction of singularity deter- 
mined by equations 13 and 15, the sin- 
gular point lies at 


s=Re? (16) 


The singularity given by equation 16 is 
then subtracted from X(s) given by 
equation 2 and the process is repeated 
until all singularities are removed. 

From example 1, Bx41=1, and 


Ee 
dy = (*) ea (17) 
n=0\ / 
Then 
ay 
— jng\1/k _ 7 — 50)|1/k 
=|1+e*| (18) 


It is seen that the right-hand member of 
equation 18 is equal to 2 for 9=0. 


Application to a Nonlinear Feedback 
System 


Consider the nonlinear feedback system 
of Fig. 3. The differential equation 
governing the dynamics of the system is 
given in the following 
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dd 
Feta +x2=r(t) 


at? ce 


The forcing function or input as shown 
Fig. 4 and given by 


r(i)=e"—2e * (20) 


Assume that the initial conditions are: 
%(0)=Oandx’(0)=1. Substituting equa- 
tion 20 in equation 19 and partitioning at 
the first linear term gives 


Ox ig —3t 
Dae ate —2e —(2x’+x+22) 


(21) 
The auxiliary equation is!? 


Dis —2e tn! ta tet) (22) 
n=0 


and the solution of equation 21 is given 
byt? 


x(t)= D7 CnOn(t) + Kot Kut (23) 
n=0 
where 
t 
Qn(t) -| rt” y(t—7)dr (24) 
0 
and 
1 ft 
y(t) = all = ds=tu(t) (25) 
27j J pp S? 


Fig. 3. Nonlinear feedback system 
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Fig. 4. Sketch of input and response func- 
tions 


where Br is the Wagner-Bromwich con- 
tour enclosing the origin and w(é) is the 
unit step function. Note that s? inside 
the integral is a special form of Z(s), which 
depends upon the form of Z(D) in equa- 
tion 3. In equation 21, Z(D)=D?; hence 
Z(s)=s*. Substituting equation 25 in 
equation 24 gives 


t - pete 
Ont Mente FD 


Applying the initial conditions gives 
Ko=0 and K,=1 in equation 23. Thus 


(26) 


pote 
x(t G fae) 
he 3 "(tin +2) * ae 
Substituting eqution 27 in equation 22 
and solving for the constants recursively! 
gives 


Q= —3, (ire (C3= —15/2, C;=31/6, Bats 


(28) 


Substituting these values in equation 27 
gives 


x(i)=t—5 P+ BR e— x8 P29) 
The 7 th term is recognized as 
(=D QE=D 
n! 
so that 
x(t) = eo) ” (30) 


The Laplace transform of equation 30 is 


=f n+1 gn il 
X(s)= oe eas (31) 


Stability Study of the Feedback 
System 


From equation 11, the modulus of the 
first singularity is 


R= lim|(—1)"*1(2"—1)|'/"=2 (32) 
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- Since the radius is not unity, application 


of equation 14 gives 
Bayi = —(—1)"(2"— 


Substituting equation 33 for Bn41 im 
equation 13 gives 


k 
d= -r(3) —1)"(2"—1)2° "7"? (34) 


It can be shown that for 6=7 


LD eg (33) 


k 
ee B\ =n i/k 
iF a7-gr—1)/F=2 (35) 
anes mG} 
Thus the singularity is at 
wr =2eo (36) 


If Bn+iR”*! is a power function s is a 


pole. 
From example 3, equation 8 gives 
1 — nee 
hGQ=a5- Se Sr (37) 


fora=2. Subtracting Xi(s) from X(s) as 
given by equation 31 yields 


¥)=X()- HL) = =) at (38) 


The singularity of equation 38 is at 
n=” (39) 


by inspection. Plotting the singularity 
hull, shown in Fig. 5, and applying the 
stability criterion, the system is stable. 

In this example, the singularities are 
poles. In simplified form, equation 31 
corresponds to 


if 
stl ey 


Taking the inverse transform of equation 
40 gives 


X(s)= (40) 


x(t) te 7! (41) 


which is the solution of equation 21 in 
closed form. The first derivative is 


x'(th=—e ‘+26 " 


At t=0, equations 41 and 42 give respec- 
tively «(0) =0 and x’(0) =1. 


(42) 


Conjugate Indicator Function and 
Diagram 


So far the tranform pairs x(1) and X(s) 
have been considered. Suppose that ¢ 
be a complex variable denoted by 


90 


t=re (43) 


where r=|tl. Consider Fig. 6. The 
growth of x(t) is determined by the in- 
dicator function 


f(0)= lim 


y ie al?) 


log | «(t)| 
; (44) 


‘on |t|, which is independent of 6. 


s-plone 


Fig. 5. Singularity hull of X(s) 


This diagram involving 0 and ¢ is calle d th 
conjugate indicator diagram. The fun 
tion f(0), which is also called the Phra; 
mén-Lindelof function, determines th 
different behaviors of x(é) in the differer 
directions of fas 7. To unders 
the meaning of f(#) consider 


x(t) Zs (at so) t 


Substituting equation 43 in equatior 
and simplifying gives 


| x(#)| = <7 (@ cos 6) sin 8) 
The behavior of |x(#)| depends first of ¢ 


as t+ © the growth of x(¢) depends ot 
sign of the factor involving 6. That 
say, along different rays of ¢ (as def 
by equation 43), the value of @ detern 
whether (a cos 6 —b sin @) is negatiy 
positive. If it is negative it is clear 
lx(#)|+0 as r>o. It is also clear 
there is a range of @ for which the fe 
(a cos 6 —b sin @) is positive. Int 
cases |x(é)| will be unbounded as 7 t 
to infinity. The function f(6) is th 
measure of whether \2(4)] is stable or 
It may be noted here that the prope 
of f(@) are the basis for the developt 
of the stability criterion given in the 


going. 


Imt- ; t-plans 


Fig. 6. Conjugate indicatee diagrar 


fice we deal with positive real time, 
nly need consider results of f(9) for 
in Fig. 6 if we are interested in the 
ility of a system but not the ‘degree 
ability. It is possible to utilize the 
‘mation obtained from f(0)40 to 
fin information about the degree of 
ility when @=0, as discussed in the 
jsection. 


prmining Degree of Stability 


ferring to equation 46, f(6) is roughly 
as follows: Since 


(2)| 


=a cos 06—6 sin 6 (47) 
for large r: 
@ cos @—6 sin 6 (48) 


seen that if the system is stable then 
=0, f(@)~a<0. An angle 6=0,+0 
s such that f(6,.)>0. We now show 
to find this @, graphically from the 


ppose that the CSH is in the LHP, 
own in Fig. 1. The critical angle 9, 
is the direction of the ray which 
2 projection of a singularity at the 
This angle 6, measures the 
tion when first the system has an un- 
ded complex time response. To 
mine this angle first plot the CSH in 
-plane. Construct the perimeter of 
SH by connecting points of singu- 
with straight-line segments. Draw 
aight-line tangent to the perimeter 
nown by b in Fig. 1) from the ori- 
Next construct a line or perpen- 
ar to line }0. The inclination 
of line Or is 6. 

peat this procedure for the tangent 
in the lower-half plane w<0. The 
so obtained is the negative critical 
, to be denoted by 6,1. The sum of 


these angles is the total critical angle 67. 
This leads to the definition 
6p 


ee 
7 


(49) 


as the degree of stability. 

This equation gives a first-order measure 
of the degree of stability. Since 6,< T, 
nS1. When »=1 the system is non- 
oscillatory and exponentially damped. 
To compare two such systems a second- 
order measure of the degree of stability 
might be the magnitude of the displace- 
ment of CSH along the o axis. This 
would be the magnitude of the distance 
of the nearest singularity to the origin 
lying on the real axis «. A third meas- 
ure would be the distribution of the pro- 
jected singular points along the critical 
ray. 


Conclusions 


In this paper we have presented some 
aspects of a stability theory and a gen- 
eralized stability criterion for the analysis 
of a certain class of nonlinear systems. 

In nonlinear systems the stability is 
markedly dependent upon the boundary 
values (the initial conditions, for ex- 
ample). It is also dependent upon the 
magnitude of the driving function. In 
general, the range of amplitudes of the 
driving function and the range of initial 
or boundary conditions would have to 
be specified. A stability criterion for non- 
linear systems must involve all these con- 
ditions in addition to the characteristics 
of the system. The proposed criterion 
based upon the concept of the convex 
singularity hull contains all these condi- 
tions automatically. 

Because of the conformal properties of 
the CSH and the conjugate indicator dia- 
gram, these diagrams can be superimposed 


and relations between the frequency and 
time domains of a given nonlinear system 
obtained immediately. 

In references 1, 2, 3 it is shown how a 
nonlinear differential equation such as 
equation 3 can be transformed to a re- 
currence equation. The latter can then 
be solved in open form because of its con- 
volution and symmetry properties. In 
many instances the recurrence solution 
can be converted to a closed form by in- 
duction. Then the present method is 
immediately applicable. It appears from 
the Taylor-Cauchy inversion transform 
reported in reference 3 that it may be 
possible to ascertain the stability of a non- 
linear system directly from the recursion 
formulas. 

If the closed form of h,, in equations 1 
and 2 is not known, it may be possible to 
obtain roughly the stability of the system 
by taking only the first few terms of the 
expansion of x(#). In this case we must 
conceive a new function which majorizes 
x(t) at least in the known truncated por- 
tion. However, caution must be taken 
to avoid serious errors if the unknown or 
discarded part is not majorized by this 
new function. 
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Discussion 


Tsang (University of Mississippi 
ity, Miss.): The authors are to be 
itulated for dealing thoroughly with 
orny subject of stability of nonlinear 
In fact the criterion is applicable 
function, not limited to those derived 
nonlinear systems. 
-concept of the CSH is introduced. 
e LHP is not a concave region, it 
that a simpler statement ‘“‘all the 
irities must be in the LHP” should 
ficient for our purpose. 
re will-be considerable difficulty both 
iving the Laplace Transform from 
nlinear equations with their boundary 
ions, and in determining the singular 
from a Laplace Transform in inverse 
series of s. Equation 41 from 
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© 


equations 19, 20 with x(0)=0, x’(0)=1 
is not typical for nonlinear systems. It 
would be very enlightening if the authors 
will apply the method in detail (first obtain 
a power series, then its Laplace Transform 
term by term and the singularities) to 
equations 19, 20 with the boundary condi- 
tions, say x(0)=—1.9, x’(0)=0.2. 


Y. H. Ku and A. A. Wolf: The authors 
would like to thank Professor N. F. Tsang 
for his interesting comments on our paper, 
and take this opportunity to thank every 
one who discussed this paper with us at 
its presentation. 2 ee 
With regard to Professor Tsang’s com- 
ment on the singularities, reasons shall be 
briefly given for the necessity of convexity 
of the singularity hull. The rate of growth 
of an entire function in different directions 
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in the complex ¢-plane can be specified as 
indicated in reference 2 of the paper by the 
Phragmén-Lindel6éf indicator function 


JO 
fi(@) = fim 22 12e — | (50) 


where p denotes the order of the function 
and ¢=reJ®, It is known from the theory 
of entire functions that if p<1/2 for x(t) 
and if x(t) is not a constant the minimum 
modulus is everywhere unbounded. That 
is to say that x(t) denotes the response of 
an unstable system. In reference 2 of the 
paper it is shown that for certain class of 
physical systems p is unity or approximately 
so. Nonlinear physical systems discussed 
in this paper are governed by functions of 
exponential type and order one. 

Putting p=1 in equation 50 gives the 
indicator function 
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— log | x(re?®)| 
a a i (51) 


fh(8) = 
eae 
for this class of physical systems. 
The indicator function is everywhere 
bounded by the type of x(t), namely r. 
The type 7 of a function x(¢) is given by 


__ log M(r) 
== iar 


mo Pe 


(52) 


according to reference 2 of the paper, where 
M(r) denotes the maximum modulus of 
x(t) and r=lel. 

Hence 


M(r)= max | «(4)| (53) 


where D, is the domain of convergence of 
x(t) and r belongs to D;. From equation 
53 it is evident that 


M(r) 2 |x(¢)| 


for all r belonging to D,;. Thus on sub- 
stituting equation 54 into equation 52 and 
comparing the result with equation 51 
one obtains the simple relation 


ONG 


(54) 


(55) 


Geometrically equation 55 signifies that 
fh(@) is associated with a convex domain. 
Because of the conformal behavior of the 
complex ¢-plane and the complex s-plane, 
the distribution of singularities is also 
convex. Hence the concept of a CSH is 
important and useful. 

This approach can be extended to those 
systems giving rise to nonconvex singularity 
hulls by suitably breaking the nonconvex 
sets into the sum of convex sets and then 
applying to each such convex set the theory 
presented in this paper. 


In answering Professor Tsang’s next 
comment, the example given in the paper 
shall be discussed with general boundary 
conditions: 


x(0)=a=Ko; «’(0)=b=Ki (56) 
where a and J are real constants. Equation 
23 in the paper gives then 
yrr2 
(57) 


10-3) & ane tet 


The general recursion formula for Cys 
given by 


Sie 26a. 
C= at —98)"] 
nN! nN 
26Cn ———— —aC, — 
(ae 
bCn1© — Cou 0 —2a 
Cn—3 
2 —a2C. (i) ees 
(n= Wea 
QabCn_1© —b*Cn_2™ (58) 
where 
Ca = 1 when m=0 (59) 
O otherwise 
Cait 
n-4 
Ci Cn Ka 


(60) 


When a=0 and b=1, one obtains the 
recurrence relation (see equation 60) 
from which equation 28 in the paper is ob- 
tained. 


- term is known as soon as the C§ 
‘specified. It is even possible to buile 


Notes on Bridged-I Complex Conjugate 


Compensation and 4-Terminal 
Network Loading 


PRAPAT CHANDAKET 


NONMEMBER AIEE 


ANY servo systems possess open- 

loop transfer functions which con- 

tain one or more pairs of complex con- 
jugate poles. Often a pair of these com- 
plex conjugate poles will lie close to the 
imaginary axis and will cause the un- 
compensated system to possess a rela- 
tively poor damping ratio or damping ex- 
ponent. Usually, satisfactory compensa- 
tion of such systems can be achieved 
through the use of tachometer feedback 
or the conventional compensators possess- 
ing real axis singularities. However, 
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ALLEN B. ROSENSTEIN 


MEMBER AIEE 


there are occasions when complex con- 
jugate singularity compensation would 
possess distinct advantages. This is par- 
ticularly true if the compensating com- 
plex singularities can be introduced by a 
simple R-C (resistance-capacitance) pas- 
sive network with zero d-c attenuation. 
Evans! has shown how the simple 
bridged-T network can be used with a-c 
servos to produce a phase lead. In this 
article, the application of the elementary 
bridged-T network is developed for d-c 
system compensation. The locus of 
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_ Control Systems Committee and app 


fee enes 
| (n-1)(n—2) 


~ Cone 


Cn—w 0 


It is always possible to solve equatior 
and 60 recursively; i.e., for ~=0, 1, 2, 
Suppose that this is done, then the 
term shall either be recognized or 
In the former case one proceeds as in 
paper. In the latter case the form 
Co, Ci, Co, ... might be examined 
if it takes a particular form. One | 
form might be 


where A and a are constants, as sug: 
by the nth term given in the pap 
this should prove to be unfruitful, i 
a majorizing or asymptotic solution 
Cn might be obtained. } 

In conclusion, it is seen that it 
straightforward matter to obtain rect 
relations which can be solved. It is ree 
nized that some difficulties may 2 ise ¢ 
getting the closed form of the mth 1 
Further efforts are being made it 
direction. It may be remarked that 
method of stability analysis can be re 
so that one can synthesize a no 
feedback control system with a pre: 
frequency response. In this case t 


the system a given degree of stabil 
given settling time, and other rel 
characteristics. 


roots of the transfer function 
bridged-T network is presented in d 
sionless form. It is found that 
transfer function of the brid 
possesses two zeros, one or the 
of which may be adjusted to fall upo: 
desired point in the left half-plat 
cluding any complex conjugate 
The bridged-T also produces a 
poles which fall upon the nega’ 
axis on either side of the real ¢ comp 
of the pair of zeros. 

A set of curves has bees dew 
which allows very rapid synthesis 
unloaded bridged-T network wil 
compensating complex conjugate 
falling at any desired left half 
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1. Common types of bridged-T con- 
figuration 


mts. Both types of simple bridged-T 
work can be synthesized from the 
ne graphs. 
fo investigate the effects of loading 
pn the bridged-T network, a general 
hod has been developed to observe 
locus of the transfer function poles 
i zeros of any four terminal passive 
work with varying load termination. 
e effects of three different types of 
dged-T loading are investigated. It is 
nd that these loads do not affect the 
ation of the zeros but can be employed 
change materially the location of the 
mpensating poles upon the negative 
axis. 
he application of the loaded bridged-T 
work to a position servomechanism, 
h as a fin control servomechanism 
ose open-loop transfer function along 
Ih the airplane dynamics function, 
itains complex poles in the left half- 
ne, is found to offer material improve- 
nts in the over-aJl system perform- 
es. 


alysis of the Bridged-T Network 


ANSFER FUNCTION 


sridged-T R-C networks fall in two 
es of configurations; See Fig. 1. The 
work in (A) will be referred to as 
lged-T type Ja, and the other as type 
_Bridged-T type Ja has been ana- 
sd by Evans! and others. Evans used 
t locus techniques to determine its 
zero configuration and found it to 
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be characterized by two poles and two 
zeros. These poles are always real and 
negative, but the zeros can be made nega- 
tive real or complex, with negative real 
parts depending on the ratio of capaci- 
tances. 

The root locus technique will be used 
to show the movement of the poles and 
zeros of the transfer functions of these net- 
works. The bridged-T type Ib is selected 
for illustration. With zero input source 
impedance and infinite output load im- 
pedance, the voltage transfer function of 
the network type Jb is found to be 


oat S. Ht +2R, Cs +R,RoC2s2 
Ein 1+(2Ri+R2)Cs+R1RoC2s? 


(1) 


Locus OF POLES AND ZEROS 


The equation for the zeros of equation 1 
is given by 


1+2RiCs+RiR2C%s?=0 (2) 
which can be arranged as 


1 +2R, Cs a 


Riki? (3) 


The pole equation can also be arranged 
as 


iL +2R;, Cs, 


RoCs(1+RiCs) @) 


If 1/R2 is assumed to be the loop 
gain and allowed to vary from zero to 
infinity, the root locus plot of equation 3 
will give the zero configuration for the 
network, while the root locus of equation 
4 will yield the pole characteristics. 
These are shown in Fig. 2. 

From Fig. 2, it may be seen that the 
zeros of the transfer function can be made 
complex or real in the left half-plane while 
the poles are always negative real when 
R> is varied. 


SIMILARITIES 


The transfer function of the network 
type Ja with zero input impedance and 
infinite output load impedance is found 
to be 


Eout 1 +-2RGs+R?CiCs* 


= 5 
Bin 1+R( Ci +2C2)s+R?2C, Qos? ( ) 


Comparing the transfer functions of 
both networks, equations 1 and 5, it is 
obvious that they have similar character- 
istics. If R, Ci, C2 (type Ia), equation 4, 
are replaced by C, R, and R; respectively, 
equation 4 becomes the transfer function 
for type Jb as shown by equation 1. A 
single graph giving the pole-zero laws for 
both types of network in normalized co- 
ordinates would be a convenient design 
tool. Such a graph is developed in the 
next section to allow the rapid location 
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(B) 


Fig. 2. Root loci of zeros and of poles of 
transfer function versus loop gain 1/R> 


A—Of zeros 
B—OF poles 


of a pair of complex conjugate zeros 
anywhere in the left half-plane. 


Normalized Characteristics 


For design convenience, a graphical 
method is developed in this section. It 
should be made clear at the beginning 
that all derivations are given with no 
loading. It is desired to design the 
bridged-T networks with specified con- 
jugate complex zeros. The poles are de- 
dependent on the zero locations. How- 
ever, loading influences the pole-zero 
location of the transfer function. This 
effect will be discussed in the next section, 
where it will be used as part of the 
synthesis procedure. 


NorMAL DESIGN CHART DEVELOPMENT 


As shown, type Ja and Jb networks are 
equivalent to each other. The same 
transfer function can be obtained just 
by switching their resistances and capaci- 
tances. Therefore, for the development 
of the design charts, only network Jb will 
be analyzed. 

As shown in Fig. 2(A), the locus of com- 
plex zeros is along the circle centered 
at —(1/2R,C). If equation 2 is re- 
written in terms of its damping ratio and 
undamped natural frequency (w,) as 
follows: 


by 
RRC RC 


=n? +2fwnsts?=0 (6A) 


1+2RiCs +R RoC2s? = oe 


then 
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Fig. 3. Design chart for bridged-T type la and 
Ib taken from normalized equation 5 
Eout _ 1+92s+as? 
Eis) 14+-@-+a)s+as" 


SCALE FACTOR: TYPE Ia: I/RCo 
Ib: VCR, 


(_) BRACKETS INDICATE RATIO OF C,/Cp OR Ro/Rj. 


(1.7) 
(1.5) 


(1.7) 


NOTE: ZEROS FALL UPON REAL AXES WHEN RATIO=# I. 
RATIOS LESS THAN | ARE NOT SHOWN. 


io 
a 


iy 


-18 -I7 -16 —-I15 -I14 —13 

undamped frequency =a = et 
CVRiR2 

damping ratio={=V. Ri/Re (6B) 

and 

wn/f=1/RiC (6C) 


The damping ratio ¢ is seen to depend 
upon the ratio, R:/R;. To normalize 
the locus of zeros, Fig. 2(A), a scale 
change can be introduced to cause the 
circle diameter to be equal to unity or, in 
other words, to cause 1/RiC=1. With 
R,C=1, equation 2 is normalized and can 
be rewritten 

1 2 Ry Ri 
RAR RG Paces sth eas 

=n? +26nz0n25 +5? =0 (7) 


from which it may be seen that 


normalized natural frequency 


j beng 
a SA 
Baa eM 


=Wnz= 


normalized damping ratio 


Sy gts hid 
bua VR /Re= (8B) 
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Ss 


10 
POLES P,P, 


30:9 3-010 uae ON ena O16 


if itis assumed that a= R;/R. 

But from equations 6(B), 8(A) and 
8(B), it is noticed that the original damp- 
ing ratio, the normalized damping ratio, 
and the normalized natural frequency 
of oscillation are all equal. 
(eee Cine (9) 

Va 

For complex conjugate zeros, it is neces- 
sary for Ry to be greater than R;. Nor- 
malizing equation 2 into equation 7 by 
setting 1/R,C=1, the locus of zeros can 
be plotted, Fig. 3, as the ratio of R:/R; 
is varied from unity to infinity. Sim- 
ilarly, equation 4 can also be normalized 
with the use of the scale factor R,C 
and the locus of poles plotted in Fig. 3 
as the ratio R:/R; is varied from unity to 
infinity. Note that the plots could have 
been presented in terms of Rj/Ro, but 
it has been considered desirable to use 
large numbers rather than small decimals. 

The ability of the ‘bridged-T network 
to synthesize any left half-plane complex 
zero can be realized from Fig. 3. Here it 
can be seen that complex zeros can be 
readily placed with any desired damping 
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—0'5 


(2) (2.5) 
(3) 


(2) 


ratio between zero and unity. Furth 
since the damping ratio is invariant wi 
scale change, the designer has only 
change scale by selecting the proper val 
of 1/RiC to adjust the complex zet 
to any desired values on the required 
line. 


works of both types Ja and Jb. 

numbers in the bracket indicate 
ratio C,/C2 for network Ja or Ro/R; | 
network Jb. Direct readings of poles a 
zeros can be obtained only under 1 
condition that 1/RC, (network Ja) 
1/CR, (network Jb) equal unity. 1 
exact location of the poles and zeros ¢ 
be obtained by multiplying the chi 
reading by 1/RC2 or 1/CR;. These f 
tors are therefore the mentioned se 
factors for the determination of the act 
pole-zero locations. On the other hai 
the desired locations of the zeros might 
given and the problem would then be ¢ 
of determining the required network. 
this case, the given value of ¢ would p: 


duce R2/R; or C\/C2 from the chart. 7 


ratio between the actual co-ordinate 
the zeros and the normalized co-ordine 


ae. 


Jury 


from normalized equation 5 


Eout _ 1+9s+as? 
Ein = 1+(2+a)s+as? 


SCALE FACTOR: TYPE Ia=1/RCo 
| TYPE Ib=1/CR, 


d produce the scale factor which 
1/RC2 or 1/R,C. 

r compensation in many practical 
ications, one must place the com- 
“zeros close to the imaginary axis. 
obtain more accuracy, an expanded 
t of the region in the vicinity of the 
inary axis has been provided in Fig. 


IGN CHART APPLICATION 


vo examples of the application of the 
ged-T charts are given: 

sample 1: To determine of the volt- 
transfer function of the network shown 
ig. 5(A). This is a type Ja network 


/ RC.= 10, and C7 Co— 10% 1/RCy 
mes the scalefactor. From the chart 
ig. 4: 


ys: —0.1470.3 

0.09, —1 11 

fultiply singularities by 1/RC, to ob- 

the correct transfer function. 

_(s+1473)(s +138) 
(s+0.9)(s-+11.1) 

xample 2: To design a bridged-T 

: Ia network to give a pair of complex 


s at ¢=0.2 and w,=2, and then find 
its voltage transfer function. 


¥ 1959 


4. Expanded chart for small damping ratio taken 


( ) BRACKETS GIVE THE RATIO OF C)/Co OR Ro/R;. 
CJ BRACKETS GIVE THE Cn FOR THE GIVEN RATIO. 


lied 0) 
POLE P5 


ala) 


From equation 6(B): If ¢=0.2, then Ro/R; 
orGy/C,=25 
From the expanded chart, Fig. 4, for 
Ci/C.=25 
Zeros: —0.04+70.195 
Poles: —0.0384, —1.042 
From equation 6(C), 
URE =a7)/C=2/02=10 
Actual pole-zero location: 
Zeros: —0.4+71.95 
Poles: —0.384, —10.42 


Transfer function 
HS +0.4+71.95)(s-++0.4 —71.95) 
~ (s-++0.384)(s-++10.42) 


Network type Ja can be obtained by 
selecting its components under the condi- 
tion that G/C,=25 and 1/ROQ=10. Tt 
one selects R=1 megohm, then C,=0.1 
microfarad and (C,=2.5 microfarads. 
This network is shown in Fig. 5(B). 


SUMMARY OF PROPERTIES OF UNLOADED 
BRIDGED-T 


From the previous analysis, the proper- 
ties of both types of networks without 
loading may be summarized as follows. 
Refer to Fig. 1(A) and 1(B): 

1. Product of poles=product of zeros (see 
equations 1 and 5) 
2. Poles are always on negative real axis. 


3. Zeros can be made negative real or com- 
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(5)C0.446] 


(50)L0.142) 
(60)C0.129.] 
(75)Lo0.116] 
(100) Lo.1 


plex according to the ratio of C\/C: (Ja) or 
R2/R; (Ib). 
Network Ia: 
C,/G:>1, complex 
C,/C,=1, double, negative real 
Ci/CQ<1, different, negative real 
Network J6: 
R2/Ri>1, complex 
R:/Ri=1, double, negative real 
Rs/Ri <1, different, negative real 


4, One pole is always closer to the origin 
than the real part of the zeros. 


5. On the complex plane, the locus of 
zeros is a circle. 


6. No d-c attenuation. 


7. For very large Ci/C, or R2/Ri, one pole 
and two zeros approach the origin; another 
pole approaches —1/RC, for type Za and 
—1/CR, for type Ib. 


APPLICATION 


The property g is sometimes very use- 
ful in stabilization of a certain type of d-c 
servos, such as those positioning types 
with an inner loop pole at the origin and 
two more poles on the negative real axis. 
Suppose the open-loop transfer function 
of this system is K/[s(1+-715)(1+795) ] 
In this case, the pole at the origin will 
associate with the two complex zeros and 
the one negative real pole which are all 
close to the origin; see Fig. 6. 
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(B) 
Fig. 5. Network configurations. A—For 
example 1. B—For example 2 


If RC2 (for type Ja) or CR, (for type Ib) 
is made very small, then the other pole 
of the bridged-T transfer functions will be 
very far from the origin. It will there- 
fore not be significant in controlling the 
system roots. Those singularities which 
group together very close to the origin will 
have little effect on the two remaining 
poles —1/7; and —1/72. Consequently 
the root locus with a break between these 
two remaining poles will be practically ver- 
tical and one can increase the gain to a very 
high level while maintaining a large damp- 
ing exponent for the predominant roots. 
The main limitation for this application 
of the bridged-T is the roots near the 
origin. However, if the gain can be 
made very high with the system stability 
still satisfied, then the residues due to the 


roots close to the origin will be very small, © 


and will contribute very little to the closed- 
loop transient response. This example, 
however, is somewhat oversimplified and 
must be used with caution. The pres- 


—— 

i REaatTYPE'L,) 2) ee 
T> na 

OR-1/CR, (TYPE I,) 


RBI 8 23 o 


pe | 
Nea 


EXAGGERATED Sot 


(A) 


Fig. 7. Transfer function G(s) of four terminal passive networks 


ence of the zeros in the vicinity of 
the origin will introduce a low-frequency 
notch in the open loop transfer function. 


Effects of Loading Impedance 


It is known that, if a 4-terminal net- 
work is terminated by any load imped- 
ance, the pole-zero configuration of its 
original unloaded transfer function will be 
affected and one will obtain a new con- 
figuration. It was desirable, therefore, 
to develop a technique to observe the 
change due to this loading effect. For- 
tunately, this technique can be used 
as part of the synthesis procedure in the 
bridged-T design. Application of this 
concept will be illustrated in the stabiliza- 
tion of a control loop containing airplane 
dynamics in the following section. 

Generally, this method can be applied 
to any 4-terminal network to observe the 
locus of the pole-zero change of the orig- 
inal transfer function when a variable 
load impedance is connected across its 
output terminals. It can be used, as well, 
for adjustment of the configuration. 


THEORY 


The method. employs Thevenin’s 
theorem in conjunction with root locus. 
The pole-zero change of the network trans- 
fer function due to loading can be obtained 
by the following steps: 


1. Having determined the original un- 
loaded system transfer function, use 
Thevenin’s theorem to find a new transfer 
function with load impedance Z(s) being 
placed across the output terminal. 


2. Convert this function into an equiv- 
alent feedback control system. 


Fig. 6. Root locus plot of type | d-c servo with bridged-T network 


at extreme conditions as its compensating device 
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G(s) 
E 


Fig. 8. Transformed circuit by Thevenin’s theorem 


[Lee 


(8) 


38. Apply root locus technique to 
equivalent control system in step 2. 
roots of the locus plot will now slow - 
changing poles and zeros under the va 
able load conditions. 
The three steps will be analyzed as fi - 
lows: Suppose 
G(s) is the original unloaded voltage no 
function. 
Z(s) is the load impedance function. 
G(s) is the new transfer function including 
load impedance Z(s). 
From Fig. 7(A) 


EX) 
Ex(s) 


From Fig. 7(B) 


Ey'(s) 
E (5) 


=G(s) 


=G"s) 


By Thevenin’s encore Fig. 7(B) ¢ 
be transformed to Fig. 8. The nm mn 
transfer function is 


Ex'(s) _ 


Gi"(s)= Eis) 


G(s)Z(s) 
Z(s)+Z(s) 


where Z'(s)=impedance looking be 
into the output of network G(s) witht 
input short-circuited. 

Suppose G(s), Z(s), Z’(s) are rewritt 
as the quotient of polynomials as follov 


(I+ LYSII+L) 14 
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40J 190] }00y = “O] “B14 
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REAL PART OF ZEROS 


Fig. 17. Root loci for 


s(1 + 2RiCs) 


i 
Re Gt iC bBs) 


Note that Z)’/(s) =G(s) by the property 
of all four terminal passive networks. 
Then, from equation 12 


G(s) 
1+2"(s)/Z(s) 
Then the poles and zeros of the new 


function G’(s) can be easily obtained as 
follows: 


G'(s)= (14) 


1. Zeros of G’(s)=combined zeros of G(s) 
and Z(s). 

2. Poles of G’(s) are the roots of [Z’(s)]/ 
[Z(s)]=—1, which can be found by the 
root locus technique. 

3. D-c attenuation=G’(0) which can be 
obtained from equation 14. 


Finally, G’(s) can be written as 


(1—s/ai)(1—s/az).. .1—s/an) 


IIA) (1—s/n)(1—s/re)...(1—s/rn) 
(15) 

where 

1, do, ...@y are zeros of G(s) and Z(s) 


Zs) _ 
[Z(s)] 


1, T2, ...?%m are roots of 


LOADING ErFrect ON Type Ia 


To illustrate the application of this 
method, three types of load impedance 
will be applied first to bridged-T type Ja. 
These are pure resistance, pure capaci- 
tance, and a combination of resistance 
and capacitance. 

Case Ia: Z(s) is pure resistance; see 
Fig. 9. Suppose 


1+as+bs? 
(1+a1s)(1+ a5) 
ZAG Sey 

Z'(s) is found to be 


2R(1+RCs/2) 
1 +R( Ci +2C, )s +R2C, Qs? 


G(s) = 


Z"s)= (16) 


or 
2R(1+RCs/2) 
Z(s)=—— 
OE rteas\(l-axs) 
Zs) 2R_(1+RCs/2) _ 
Zs) Re (1+eais)(1+ass) 
A sketch of the root loci for poles of 
G’(s) is shown in Fig. 10 by equating 


—1- (17) 
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Z'(s)/Z(s)=—1. 
Finally, 
R, 1+as+bs? 
2R+R, (1 —s/n)(1 —s/r2) 
Analysis of Bridged-T type Ja with 
load resistance R;, shows the following 
properties: 
1. D-c attenuation=R;/(2R+Rz2). 
2. Zeros are the same as original. 


3. Both poles of G’(s) are real and these 
move to the left if R; is decreased. 


(18) 


G"(s)= 


4. Limit of poles movement for the one near 
the origin is —2/RC,; or twice the real part 
of zeros. The other pole will move towards 
infinity along the negative real axis. 

It is of interest to note that, with the 
introduction of d-c attenuation, the poles 
of the original bridged-T transfer func- 
tion can be pulled farther into the left 
half-plane with the pole at —1/a, and 
are even capable of moving to the left of 
the real part of the zeros. 

Casella: Z(s) is pure capacitance; see 
Fig. 11. Assume that G(s) and also 
Z'(s) are the same as in case Ja but that 


Z(s)=1/C,s. Then 
Zs) _ s(1+RCs/2) oy: 
Z(s) Reger ie aise 


The root loci of Z’(s)/Z(s)=—1 gives 
the locations of poles in G’(s). The loca- 
tions of the new poles of G‘(s) depend on 
C1 if G(s) is fixed. The pole movement 
versus the loop gain C;is shown in Fig. 12. 
G'(0)=1 
Then 
1+as+bs? 
G'(s)=———* _. 

OG sin —s/) 

The analysis of bridged-T type Ja with 
pure capacitive load C; shows the follow- 
ing properties: 


(20) 


1. It contains no d-c attenuation. 
2. Zeros are the same as the original. 


3. Both poles of G(s) are real and move to 
the right if C; is increased. 


4. When C; approaches infinity, one pole 
will approach the origin and the other pole 
approaches —2/RC,, which is twice the real 


part of zero. 
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x O 
I [ I 


— 3A 
Fig. 18. Root loci for 


s(1 +2RiCs) 
(1 +1s)(1 +608) + RCs) 


Case IIIa: Z(s) is a series combina} 
tion of resistance and capacitance; 
see Fig. 13. G(s) and Z’(s) are taken the 
same as the previous two cases. 
1 +RiCis ‘ 

Cis x 
Zeros of G'(s)=(1+as+bs?)(1+R:C,S) 
Poles of G'(s)=roots of equation 

s(1+RCs/2) 3 
(1+ ais)(1+e2s)(1+RiCis) 


Zs) = 


2RCi 


Root loci for expression 19 are shoy 
in Fig. 14. ; 

If R, and C, are arbitrarily selected, a 
zero in G’(s) will be fixed at —1/Ri 
-also, the poles of G’(s) are fixed. 
viously, this network has no d-c attent 
tion. Suppose the roots of the syste 
are fixed at 1, 72, and r3; then 


(1+as+bs?)(1+RiCis) 
(1—s/n)(1—s/re)(1—s/rs) 


To make this analysis applicable, 
should be emphasized that one negatt 
real zero may be added at any point ¢ 
pending on the selection of R; and 
Poles can be adjusted by varying C; wi 
the condition that R:C; be kept 
changed. Under this condition C, (g 
equation 21) can be considered as t 
open-loop gain of the system. ) 

Analysis of G’(s) of bridged-T type 
with series resistance R; and capacitan 
C, as load-impedance shows the followi 
properties: : 


G 


G'(s)= 


1, There is no d-c attenuation | 


2. Contains three poles and three ze 
(two original and one negative real zero 
—1/RiC,). 


3. For a specified position of new ze 
poles can be adjusted by varying C, w 
R,C, maintained constant. 


LOADING EFFECTS ON TypE Ib 
Analysis of bridged-T type Ib w 
load impedance can be obtained by 
same technique used for type Ja. Int 
section, only the results and useful 
pressions will be given. Since the use 
these networks will usually concern ec 
plex zeros, the zeros of G(s) discuss 


Jury 1 


igate complex; see Fig. 15. 


a l+a’s+b’s2 
(1+fis)(1+2s) 


All 8;, B:, a’, and b’ are fixed and will 
F used for all cases in later discussion in 


Ss) (23) 


G(s) =1/R,C. 
a R(1+2R,Cs) 
(1+8:s)(1+/25) 
Case Ib: Z(s) is pure resistance (R,). 
pros of G’(s)=zeros of G(s), and poles 
G’(s)=roots of 
(1+2R:Cs) | 
(1+6is)\(1+é:s) 
See Fig. 16 for root loci of equation 25. 
R, is fixed and the roots are at 7, and rs, 


(s) (24) 


(25) 


eRe (1+a’s+b’s?) 
Re+R; (1—s/n)(1—s/re) 
Case IIb: Z(s) is pure capacitance C,. 
tos of G’(s)=zeros of G(s), and poles 
G’(s)=roots of 
. s(1+2R,Cs) 
a — || 27 
*(1+6i5)(1 +s) vil 
See Fig. 17 for root loci of equation 27. 


C, is fixed, the roots are at 7; and 7, 
en 


(26) 


__ (1 +a's-+b’s?) 
Se Gs /r(d—s/n) det 


Case IIIb: Z(s) is the series combina- 
mn of resistance (R,) and capacitance 
i). Zeros of G’(s)=(1+as+bs?)(1+ 
Cis), and poles of G’(s) =roots of 

s(1+2R:Cs) : 
~~ (A+6is\1+6s\(1+RiCis) 


See Fig. 18 for root loci of equation 


—1 (29) 


K (I+10S) 
S(I+0.2S$+ 0.25 $4 


Fig. 19. Control sys- 


C(s) 
: tem for example 1 


29. If R,and C, fix the roots at 7, 72, and 
r3, then 


(a +a’s +b’s?)(1 +RiCis) 
Cl=s/rl —s/r2)(1—s/rs) 


The properties of loading effect in 
these three cases for network type Jb 
are generally similar to that of type Ja. 
The differences are only the d-c attenua- 
tion terms and the limits which the roots 
of [Z2’(s) ]/(Zs)]=—1 will approach when 
the loop gain is increased towards infinity. 
Even in case JIIb, after another zero has 
been added, the poles must be adjusted 
(as in case JJJa) by varying C; with RC, 
maintained constant. 


G"(s)= 


(30) 


USE OF LOADING EFFECTS 


The behavior of the voltage transfer 
functions of both types of bridged-T net- 
works under the three loading conditions 
have been discussed in the last two sec- 
tions. One should note that, under the 
same loading condition, the two types of 
network behaveinasimilar manner. The 
only slight difference is the position of 
one zero in the equivalent open-loop 
transfer function Z’(s)/Z(s) for the deter- 
mination of poles in G’(s). This zero 
will give a different termination of one 
root. 

For example, in case Ja, this zero is at 
—2/RC; (see equation 17 and Fig. 10) 
and becomes the terminating zero for the 
root locus from the pole —1/a;. This is 
twice the real part of complex zeros of 
G(s). In case Jb, (see equation 23 and 
Fig. 15) this zero is at the center of the 
circle passing through complex zeros of 
G(s) and the origin. Consequently, even 
though G(s) for networks Ja and Jb are 
the same, one will not obtain the same 
pole locations in G’(s) for equivalent load 
resistance R;. 


eS 21. Control system with bridged-T as compensating network G(s) 
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Fig. 20. Sketch of root loci for 


K(1+10s) 
s(1+0.9s+0.25s?) 


It has been shown that the damping 
ratio ¢ of the zeros in G(s) is controlled by 
R,/R, for type Ib and C,/C. for type Ja. 
If the G(s) in both types are the same, 
then with very small ¢ the limit to which 
the pole —1/8; (type Jb) can move away 
to catch a zero at —1/(2RC) is greater 
than that to which the pole —1/a, (type 
Ia) can move to catch a zero at —2/RC. 
Thus, if ¢ is large, then type Ja would be 
preferred. 


__-19.42 


-l2 -10 -8 -6 
=.0983 
(NOT 
SHOWN) 
ORIGINAL POLE a4 
CANCELLED BY 
ZERO OF THE 
BRIDGED-T 
-6 
-8 
10 
Ip -12 
“14 
‘ a ~ —e | 
' lie 
Fig. 22. Root loci of 
40K(s+0.1) 
s(s+0.304)(s+10.42) 
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COMPENSATING SERVO CONTROL 
NETWORK 


With a knowledge of the root locus 
technique and the results of the analysis 
given in the last two sections, the proper 
type of bridged-T can be selected to meet 
system requirements. Generally, the 
effect of loading can be used to adjust 
the pole positions of the transfer function 
of the bridged-T. However, while it is 
impossible to fix both poles at will, one 
can fix the one pole which is most signif- 
icant. The capacitive load for bridged- 
T type Jb is used for pole adjustment 
in the numerical example 2 in the Ap- 
pendix. 

In actual applications, the bridged-T 
without load should be designed first by 
consultation of the design charts. Zeros of 
this network are fixed at the desired posi- 
tions; then the complete pole-zero con- 
figuration of G(s) is obtained. Modifica- 
tion of G(s) by loading effects is made by 
selecting the proper type of loading. The 
exact pole-zero location of G’(s) is ob- 
tained by applying the root locus tech- 
nique to the equivalent control system 


Z'(s)/Z(s). 
APPLICATION 


With the developed design chart, one 
can easily place a pair of conjugate com- 
plex zeros at any point in the complex 
plane. However, their associated un- 
loaded poles cannot be controlled. For 
optimum compensation it is often de- 
sirable to move the bridged-T poles. 
This can be done by connecting the proper 
load impedance at the output terminals 
and using root locus techniques to deter- 
mine the new pole-zero configuration. 
With this method, only one pole can be 
fixed in a specified position. In many 
cases, however, there will be one pole that 
is more significant than the other, and 
this technique may still be used. 

Probably the most common applica- 
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Fig. 24. Sketch of root loci of equation 35 
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K(I+10S) cls) Fig. 23. Block 
diagram of the 
S(i+2S+.25S ) system for went 
AIRPLANE DYN. ample 9 


tion of bridged-T networks to d-c servo 
systems would be in those systems in 
which the transfer functions contain com- 
plex poles. The simplest way to com- 
pensate these types would be use of the 
cancellation technique of putting zeros on 
top of the complex poles. The system to 
be employed for illustration is the air- 
plane auto pilot. Two examples are 
given in the Appendix. In example 1, 
only the bridged-T without loading effect 
will be used, but in example 2, the load- 
ing effect will be applied as well. 

It should be stated that these examples 
are mainly intended to illustrate tech- 
niques for improving the system sta- 
bility. If all specifications of the system 
are not entirely met, additional measures 
would, of course, be used. It should also 
be recognized that in actual systems, the 
exact pole-zero locations may not be 
exactly known and will often vary with 
change in environmental conditions. 
Exact cancellation of terms is, therefore, 
seldom if ever achieved. 


Fig. 25. Root loci of the system compensated by bridged-T without idea 
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Appendix 


Numerical Example 1 


The control system shown in Fig. 19 
selected to demonstrate the use of a br 
T network. G(s) is the compensati 
bridged-T network. The open-loop cont 
transfer function, in this case, represents 
short-period longitudinal mode of an a 
plane. The fin control servo is ass 
perfect to avoid its effects on system | 
formance. Without compensation, or ( 
=1, the root loci of this system are shi 
in Fig. 20. Without compensation, © e 
cause of the small damping ratio, the sys 
tem transient overshoot will be high. j 


With Bridged-T as Compensating 
Network 


Pole cancellation will be used to improy 
the stability of the system. The « 
plex poles in the original transfer func! 
have a damping ratio of 0.2 and a nat 
frequency of 2 radians per second. 
applying the cancellation method, 
bridged-T network is designed to give cor 
plex zeros at ¢=0.2 and w,=2. Genera 
either type of bridged-T can be used. 
this particular case, the bridged-T type 
has already been synthesized in the previ 
example 2 and will be used as the comp 
sating network G(s). The complete di 
gram of this control system is shown 
Fig. 21. From the previous result . 
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e Fi(s) is the new open-loop transfer 
tion, or 


40K(s-+-0.1) 


~ s(s-+0.384)(s+10.42) 2) 


t loci of this function are shown in Fig. 


this particular case, with a bridged-T 
ompensating network, the new root loci 
more satisfactory results. The loci 
nearly vertical, as the zero at —0.1 and 
pole at the origin are very close together 
compared with the other singulari- 


= —5.244712.1 
—0.0983 
ith r(t)=unit step input, residue at 


0.537 |156° 
0.537 | —156° 
—0.0172 


=1—0.0172¢-0'0983' 44 974-4! 
cos (12#+156°) (34) 


Vith ¢=0.4, the overshoot from 7; and rz 
%*. In this case, the residue of the 
tz (with unit step input) is nega- 
and very small, so the overshoot will 
lightly less than 25%. Actually, due 
he very long time constant of rz, it will 
a long-tail effect to the closed loop 
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Fig. 26. Sketch of root loci 
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10.42 =5 —.384 

performance. However, the residue at this 
root is so small that, in practice, one can 

= _(1+0.25-+0.255?) +0.2s +0.255?) (31) neglect it. 
(1+2.7s)(1+0.096s) It should be pointed out that this analysis 
are is made on the basis of exact cancellation 
e: (1+10s) (32) which, of course, is impossible. The tech- 
s(1+2.7s)(s+10.42) nique can still be applied for if the zeros are 


still near the poles, the residues at those 
roots will be very small and they wiil not 
effect the closed-loop performance signifi- 
cantly. 


Numerical Example 2 


The application of bridged-T with load- 
ing effect will be illustrated in this example 
to improve the system performance. The 
same problem used in example 1 is selected 
for illustration but, in addition, another 
servo control unit is placed in cascade with 
the original. This system is shown in 
Fig. 23. 

With G(s)=1 


K(1+10s) 


~ 5(1+0.25-+0.25s2) x 
(1++0.166s +0.028s?) 


Fi(s) (35) 


The root loci of this system are shown in 
Fig. 24. Without compensation, it is seen 
that the predominating roots of the system 
are about to move to the right half-plane, 
even at the very low gain. It is desired to 
design a compensating network G(s) to im- 
prove the system stability. For this par- 
ticular problem, the bridged-T networks 
again seem to to be very applicable. 


Bridged-T Applied Without Load 


If the bridged-T network, which has been 
used in the last problem, is applied to this 
system as a compensating network G(s) to 
cancel the complex poles of the airplane 
function, the system stability is improved. 
See Fig. 25 for the system root loci plot. 

If the damping ratio of 0.4 is selected for 
gain adjustment, the roots of the closed- 
loop system are 


T1> i —1.75+74.15 


r3= —0.085 
a= —3.7 
= —9.39 


The open-loop gain K equals 0.446. The 
most undesirable features are the large resi- 
due of r3, which is close to the origin and the 
low gain. When a unit step input is ap- 
plied, the transient term caused by this 
root is —0.162¢~0-085", 

Due to this large residue and the very long 
time constant of about 11.8 seconds, the 
closed-loop response to a unit step input 
will have a long tail, and the system per- 
formance is not considered satisfactory. 
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Loading Effect Applied for Poles 
Adjustment 


The root loci plot in Fig. 25 shows that 
the closed-loop response can be improved 
by moving both poles at —0.384 and —10.42 
towards the origin until the loci a and a’ 
are about to change to the left (or 6 and b’ 
are about to change to the right). The 
required location of the new poles can be 
roughly estimated by checking the angle of 
departure from the complex poles. If all 
singularities are held constant, the pole at 
—10.42 has to be moved to —6 to satisfy 
this condition. If a capacitance is used 
for loading the original bridged-T, both 
poles at —10.42 and —0.384 can be moved 
to the right. 

Another improvement that may be con- 
sidered at the same time is to try to cancel 
the zero at —0.1 by moving the pole from 
—0.384. If an approximate cancellation 
at this point can be made, the calculation is 
simplified and the effect of the long time 
constant transient term will be eliminated. 

For this particular problem, the loading 
effect of case [Jb seems very suitable. The 
bridged-T type Jb used with load capacitance 
will have terminating loci from both poles 
at —5 and the origin. If the case JJa is 
used, the loci will terminate at —0.8 and 
origin. In case IJ, there is a possibility of 
obtaining both zero cancellation and another 
pole moving to the vicinity of —6, while, in 
case Ja, the pole at —0.384 will change only 
slightly. The previous results from the anal- 
ysis of loading case Jb will be applied for 
the adjustment of the poles of bridged-T 
type Ib. 

First, the bridged-T type Jb (without 
load) containing zeros at the same location of 
complex poles in the airplane transfer func- 
tion is designed: {nz=0.2; 1/R,:C=10; 
R2/Ri =25. 

If C=1yF is selected, then R;=0.1 meg- 
ohm and R2,=2.5 megohms. The pole- 
zero configuration of the voltage transfer 
function is still the same as for the net- 
work Ja which was used in the last example. 


ANGLE OF ee 
DEPARTURE 5 
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Fig. 28. Root loci of the system compensated 
by bridged-T type Ib with capacitive load 
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Next, using the capacitive load C; and 
applying equation 27 to obtain the rela- 
tion for pole adjustment yields 


2C;s(s +5) a 
(s+10.42)(s-+0.384) 


See Fig. 26 for the root loci of pole move- 
ment. 

If the root 7 is fixed at —0.1, then C;= 
2.95, and r,=—5.75. The network which 
is to be used for G(s) is shown in Fig. 27. 
The root loci plot of the system is given in 
Fig. 28. With the damping ratio=0.4, 
the results are 71, 72= —144798.25; 13, re= 
—4,5+73.55; K=2.04. 

If r(¢) =unit step input 


c(#)=14+1.576€ **' cos (8.25t+96°)+ 
0.89¢ 3°" cos (3.55t+200°) (37) 


(36) 


The overshoot in this case is roughly 
estimated at about 26%. 
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Discussion 


Theodore A. Savo (Hughes Aircraft Com- 
pany, Culver City, Calif.): The authors 
are to be complimented on their interesting 
paper, which helps to bridge the gap 
between what might be called network 
synthesis and network design. The root- 
locus method, used in the graphical study 
of the poles and zeros of the bridged-T 
network in this paper, can of course be used 
for other networks the extension being 
straightforward. If the decision has been 
made to compensate in the ‘open-loop’ 
sense and the characteristics of the bridged- 
T are satisfactory, then the content of this 
paper provides useful information for the 
rapid design of the network with and with- 
out loading. An interesting feature is that 
the design is affected in terms of the poles 
and zeros of the voltage transfer function for 
the network, and hence allows the network 
’ design to be combined with a root-locus 
study of the control system characteris- 
tics. 

The purpose of this discussion is to pro- 
vide some parallel opinions on the subject 
matter, to offer some extensions to the con- 
tent, and to indicate cases where the method 
has proved useful. The authors have ade- 
quately studied the characteristics of the 
bridged-T (types Ja and Ib) from a practical- 
design viewpoint. One additional con- 
sideration is the problem of parameter varia- 
tions. Table I(A) and (B) conveniently 
summarizes the effects of component changes 
on the voltage transfer function E,/Ej(s). 

It is interesting to note that the study 
of the effects of a load on the transfer char- 
acteristics of the bridged-T can be conyen- 
iently expressed in terms of the short- 
circuit admittance parameters (y12, yoo). 
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Table |. Summary of Network Characteristics for Parameter Variations 
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Table Il. Summary of Loaded Types la and Ib Networks 
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The voltage transfer function for the 
network with no load is 


YEA GS) (38) 


yoo 3 


(Since it would add little to the discussion, 
the minus sign normally associated with 
equation 38 has been arbitrarily omitted. ) 
The expression for the voltage transfer 
function if the network is loaded with an 
admittance Y; can be expressed as 


Ey 
—(S)= 
Ey ) 


E, “i Vie 
Ey yoo t+ Yr 


Slight algebraic manipulation of equation 
39 yields 


4 GS) (39) 


vay) ya / Yoo 
G(S)== (=F 7,.\ 1 a\ 
Ex Vir 2) |(2) 
1 ey ae aS 
+| (2 yoo 
Sy ae (40) 
TG 
14({— Jers) 


Equation 40 defines an equivalent feed- 
back control system in which the closed 
loop characteristics represent the transfer 
function of the network with load; see Fig. 
29. The open-loop gain from Fig. 29 and 
equation 40 is 


Ws 
Loop gain = — (=) 
y: 


22 


(41) 


Equation 41 can be correlated with the 
method used in the paper, namely 


Yn g2"S) 
Ye2 , Z(S) 

The characteristics of the loaded net- 
work G’(s) can now be studied by the root- 


locus method, the characteristic equation 
being 


(42) 


+ (S)=0 (43) 


where y22=short-circuit admittance func- 
tion for the unloaded network G(s); and 
Y;,=admittance of any arbitrary load. 
Specifically, the short-circuit admittance 
functions for the two types of bridged-T 
networks are, for type Ja: 


_14R(C.42C:)S+ CGR? 


Yoo (44) 
2R(1 +o s) 


and for type Jb: 


_14C(Re+2R)S+RRC2S? 
iat Rx(1+2R,CS) 
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(45) 


10,10 


8.80 


7.50 


6.20 


P; 


4.90 


3.60 


2.30 


Chandaket, Rosenstein—Bridged-T Complex Conjugate Compensation 
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Fig. 30. Type Ib bridged-T pole positions with load Y.=(RuCi+1)/Ri 
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Fig. 31. Type Ib bridged-T pole positions with load YL =1/Rt 
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Fig. 32. Bridged-T pole positions, unloaded 


ae authors have studied the bridged-T 
acteristics with three types of simple 
s. One additional type of simple load, 
onsidered, is the case of a parallel com- 
jon of one resistor and capacitor. 
his case the load admittance is 


RrCyS+1 

———————— 6 
Ry (46) 

function Yz/y.. becomes, for type Ja: 
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e zeros ‘of the transfer function with 
arallel Rx, Cz load, G’(S), are, of course, 
Same as the zeros of the network 
fer function without load G(.S). Both 
orks, however, have a d-c attenuation. 
bvenient summary of the loading 
s on the pole locations is given in Table 
) and (B) for the three cases considered 
= paper and for the additional case just 
nted. The information shown in Table 

useful in determining qualitatively 
ffects of loading (intentional or unin- 


By Bridged-T 
rk with source 


loading i! 
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tentional); however, in a specific problem, 
if the pole locations are to be determined 
exactly, a design chart of the form shown 
in Fig. 30 is appropriate for engineers in- 
terested in rapid laboratory-type design. 
Fig. 30 depicts the type 7d network with the 
parallel, Rx, Cx type of load and the nor- 
malized time constant of the load set at 
Txr=0.1 second. Fig. 31 illustrates the 
pole positions for the type 7b network with 
simple resistive loading R;. Fig. 32 gives 
the pole positions for the unloaded bridged- 
T. It should be noted that the pole posi- 
tions shown in Fig. 30 and 31 are computed 
with R,C normalized to unity; hence a 1 yh 
RiC scale factor is necessary. Similar 
comments apply to Fig. 32. 

It is logical to consider the character- 
istics of the bridged-T network for the case 
of source loading as illustrated in Fig. 33. 

The transfer function for the case of 
source loading is, for type Ja: 


Ss 1+2ROS+C,OR2S? 
GOS) = 49 
= TTIR(G42G)SER.G1S+ 
RC,C(R+2Rs)S? 
and for type Jb: 
oe {HOR CSR RsC2S? 
( MAE COR ER Lanse (50) 


R(Ri+Rs)C2S? 


As might be expected for resistive source 
loading, equations 49 and 50 reveal that the 
zeros of the network are not altered. Study 
of the pole migration using the root-locus 
method is best effected in the case of source 
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Fig. 34. Bridged-T pole locations with 
resistive source load Rs 


loading by using the open-circuit imped- 
ance functions z(S) and zu(S). The volt- 
age transfer function of the network with 
no load is 


Eo Z2(.S) 
xy Zu(S) 


(51) 


If an arbitrary source load Z, is added, then 
the voltage transfer function becomes 


212 


ZuitZs 


= (S)= (52) 


Manipulation of equation 52, in a 
manner analogous to equation 40, yields 


= (5) =e 
1 nme 


211 


(53) 


Equation 53 indicates that the poles of 
the transfer characteristic can be studied 
with a root-locus plot of Zs/z in a manner 
analogous to Yz/+o. for load impedances: 
see equation 43. 

The open-circuit driving point imped- 
ance, 2y(.S), for the bridged-T network is, 
for type Ia: 


_1+R(C420)S+OGRS? 


54 
ie GSC ORGS) 7) 
and for type Ib: 
1+C(2Ri+R2)S+RR.C2S? 
£3 = (55) 
ee 
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The resultant ratio Zs/zu(S) is (Zs=Rs) 
for type Ja: 


Zs RsQ,S(1+2RCS) 
ELS) See eee 56 
ee (1+a1S)(1+a2S) ( ) 
and for type Ib: 
Cc 
Z (5)= arscs( s) (57) 
(1+a1S)(1+a2S) 


The parameters 1/a; and 1/az in equa- 
tions 56 and 57 characterize the pole loca- 
tions for R;=0. 

Figure 34 depicts the pole migration de- 
rived from equations 56 and 57. With 
the frequency normalization used for the 
unloaded bridged-T, the gain calibration 
of the loci of Fig. 34 is quite simple, for ex- 
ample, for type Ja: 


RG=1 
C/O, = 
R 
RG -1( 2) (58) 
and for type Ib: 
RiC=1 
ies 
Ri~ 
R 


162 Chandaket, Rosenstein—Bridged-T Complex Conjugate Compensation 


- 0.9 "lo. 


Fig. 35. Design chart for bridged-T type la 
and Ib taken from normalized equation 5 
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Consequently, for a given x value which 
fixes the complex zeros, the variation in the 
pole locations can be studied by varying the 
ratios (Rs/R) and Rs/R2) respectively. In 
conclusion, it is obvious that if the net- 
work level can be designed to yield R>>Rs 
or R2>>Rs, then the unloaded design charts 
can be used. In cases where the source 
resistance is not negligible or drifts sub- 
stantially, the results shown in Fig. 34 
are useful. 

The case of both input and output load- 
ing has also been considered by the afore- 
mentioned method; however, the number 
of variable parameters becomes large unless 
special symmetry is assumed. 

In conclusion it is worth mentioning that 
the root-locus method of designing the 
bridged-T has been successfully used in 
effecting the compensation of both electric 
and hydraulic type of feedback systems. 
Cancellation-type compensation using the 
bridged-T network is particularly use- 
ful in laboratory-type design and testing 
when the characteristics of the “fixed” part 
of the system are not known exactly or 
when they vary because of nonlinearities. 
It should also be mentioned that a simple 
lead network in tandem with the bridged-T 
can be used to effectively move the low- 


frequency pole of the unloaded brid 
to a higher frequency. 


Propat Chandaket and Allen B. Rosen 
The authors wish to thank Dr. Sa 
his interesting discussion and his add 
to the work of this paper. The t 
short-circuit admittance paramete 
study the loading effect is more conve 
for the circuit designers. Dr. Savo’s 
of pole positions for the case of brid 
under no load and under loaded coné¢ 
of pure resistance and parallel R-C 
particularly useful for the designer. 
Since the bridged-T is frequently 
ployed as an a-c filter, the authors 
like to extend the utility of the meth 
veloped a little further by adding | 1e 
ratio and notch width characteristics 
unloaded bridged-T of both types. — 
These properties can be obtained 
easily from the normalized equation 1 
32, 


Eout __ 1+2s+as? { +2) ( 


Bi, 1Efoyeene (142 (2 
1 


* 


where a=ratio of Ci/C2 (type Ia) ot 
(type Id). 

We note from the equation that th 
‘uct of the zeros equals the produc 
poles, Ley ng= pipz= 1/a. 


ff the zeros are complex, then their 
malized natural frequency is wy1 /Va; 
equation 8(A). This, of course, is the 
uency in radians locating the center (in 
ig scale) of the notch, ie., the double 
ak point of the zeros. 

Phe d-c gain of the bridged-T is unity. 
Fhe notch ratio is defined as the ratio of 
a at the notch frequency to the gain at 
> frequency, replacing s with j1/Va) 
following is obtained from equation 1: 


(61) 


P= 


This shows that the notch ratio depends 
only on the ratio of C,/C; or Re/ Ri. 

For design convenience, the design chart 
shown in Fig. 3 has been redrawn with the 
notch ratios for the various values of 
Ci/C or Re/R, added in brackets; see Fig. 
35 


The notch width frequency is defined as 
the frequency difference between the two 
poles of the bridged-T. 

From equation 60 


CER hae 
2a 2a 


2a 2a 
This gives 


4 2 
Notualizel ioien Maio (62) 
a 


To obtain the real value of notch width 
frequency, one must multiply the result 
by its scale factor. Since pif2=1/a and 
the zeros fall at 1/./a, the pole-zero con- 
figuration on a log scale such as in a Bode 
diagram is symmetrical. 
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UTOMATIC FEEDBACK-control 
system synthesis occupies a rather 
yue place among the various fields of 
neering due to the indirectness of its 
roach. An antenna will be designed 
a certain radiation pattern because 
radiation pattern is called for in 
intended application. On the other 
, a servomechanism which is des- 
d to handle step inputs exclusively 
r be synthesized by the techniques of 
quist, Nichols, or Bode. This means 
the system is actually designed for a 
ency response: to exhibit certain 
nance characteristics when subjected 
steady-state sinusoidal inputs of 
ying frequency. Yet the system will 
be subjected to this latter type of 
it. The justification of this approach 
hat the experience of the designer 
yides the basis for expecting that the 
will perform properly in its 
al input environment if it is built 
a specific frequency response. Such 
nples of indirectness could be enu- 
ated almost endlessly. 

1 brief, a common characteristic of 
design techniques (including here 
design leading to writing specifica- 
s) seems to be that they focus at one 
ut a few restricted aspects of the 
em operation which by themselves 
mot necessarily significant in the 
al operating environment of the 
em. From the study of these aspects 
designer extrapolates mentally the 
2m performance to be expected 
in the actual operating environ- 
t. Guiding him in this extrapolation 
lis experience, engineering judgment, 
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Probabilistic Error as a Measure of 
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and some semiquantitative relations, the 
actual quantitative relationships fre- 
quently being too laborious to handle. 
In addition, the designer is also averaging 
mentally, on the basis of more or less 
conscious probability considerations re- 
garding the utilization of the output 
with various inputs. In other words, 
while the synthesis is immediately con- 
cerned with various special facets of 
performance under largely hypothetical 
conditions the real concern is still the 
performance in the actual environment 
during the utilization of the output. 
Thus, looking behind all the indirectness 
of automatic feedback-control design, it 
develops that, although largely sub- 
conscious, actual system design is based 
on the one ultimate measure of per- 
formance: How much error is there at 
the times when the system output is 
utilized? 

In this paper an attempt is made to 
give this measure of performance mathe- 
matical formulation which can be the 
basis of evaluating existing systems or the 
criterion for designing optimum systems. 
This measure is believed to be more 
general and to have more clear-cut 
physical meaning than existing criteria. 
In contrast to most standard techniques, 
which concentrate on isolated phases of 
performance such as transient or steady 
state, the proposed measure unites all 
of these, favoring none, in a single 
concept. 

This paper is intended as the first of a 
group of papers. Accordingly, the basic- 
philosophies are expounded here, the 
mathematical formulation and definition 
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of the measure is introduced, and it is 
demonstrated that these definitions are 
amenable to mathematical analysis. 
General solutions rather than working 
equations are presented. 


Proposed Measure of 
Performance Evaluation 


An automatic feedback-control system 
exists to convert its input into a desired 
output which may or may not be identical 
to the input. The difference of desired 
and actual outputs is the error of the 
control system. Consequently the per- 
formance of the control system is charac- 
terized by the amount of the error at 
such times as the output is utilized. 
Clearly, this verbal statement of the 
evaluation of performance contains two 
elements: 1. the amount or size of the 
error; and 2. the times when the output 
is utilized. Both of these elements are 
subject, in their nature and relative 
significance, to the environmental condi- 
tions of the specific control system. 
For instance, in some control systems 
enlarged errors are increasingly, and 
possibly cumulatively, objectionable. 
Frequency control of some power systems 
might be an example. In other cases, 
such as in an automatic milling machine, 
errors up to a certain tolerance are not 
objectionable but all errors exceeding 
this tolerance are equally and totally 
objectionable. For most control systems 
(even for the examples mentioned) the 
lines are not so clearly drawn and the 
situation is intermediate between the 
two extremes. The point made is that 
the specific environment and function of 
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each control system will set a specific 
penalty valuation on errors of a given 
size. This penalty valuation can be 
expressed for the given application in the 
form of a penalty function F(e), a 
single-valued function of the error e. 
In a more general sense, the penalty 
valuation may vary with extraneous 
parameters such as rate of flow in a 
temperature control system or with time, 
Fe, t, v1, 02...Ur), where through 2, 
are the parameters. Of course error, @, 
itself certainly is a function of time, ¢. 

The second element of performance 
evaluation, the time when the output is 
utilized, shows equally wide individual 
variation depending on the environment 
of the control system. Some control 
systems, such as a frequency or voltage 
control on a power system, operate at 
all times on essentially steady reference 
inputs and a stationary random noise 
input, the varying Joad. All times for 
utilizing the output are equally probable 
here. The other extreme would be a 
control system which makes one adjust- 
ment upon activation by a step input 
followed by one instantaneous measure- 
ment after a fixed-time interval. Here 
the time of utilization is one instant 
after a fixed-time interval which has 
probability one while all other time 
intervals have probability zero. If the 
measurement itself takes a certain time, 
all times within this measuring time are 
equally probable. In a missile, the 
output is utilized at one single instant, 
but over a number of launchings the 
time intervals elapsing from launching to 
explosion will form a probability dis- 
tribution centering around the nominal 
range, 

Essentially, then, for any automatic 
feedback-control system, the times 
elapsing from activating the system to 
all times of utilization of its output can 
be arranged into a probability distribu- 
tion, P(t), over many repeated activations 
by identical inputs of the same system or 
identical systems. 


Mathematical Formulation of the 
Measure of Performance 


The measure of performance described 
in the preceding section can be mathe- 
matically formulated as 


r= [Oem . . Ur) p(t)dt (1A) 
0 


where the symbol s is the end sigma, the 
form of letter sigma used at the end of 
Greek words. 

For one specific type of input, a 
“deterministic” input, equation 1(A) 
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clearly gives the statistically correct 
average value of the penalized error at 
all times when the output is utilized. 

If the environment of the control 
system includes a number of different 
types of inputs, possibly an infinite 
number or ensemble of different types of 
inputs, an averaging over, this ensemble 
will also be necessary to obtain the mean 
penalized error. Then the measure of 
performance takes the form 


=) F(e(t),t,01. . .Ur)p(t)dt (1B) 
0 


where the straight line denotes the 
averaging over the ensemble, so that 
equation 1(B) describes the average 
penalized error. 

Frequently the probability distribution 
of utilization p(t) will be independent of 
the type of the input. Then equation 
1(B) reduces to 
=|, F(e(t),t,v1. . .v;)p(t)dt (1C) 

The value of the proposed measure of 
performance as defined in equations 
1(A)-1(C) is conceptually the average 
value of penalized error at such times 
when the output is utilized, and for pur- 
poses of identification it will be referred 
to by either of the names “probabilistic 
error’ or ‘end sigma error.” As a 
performance measure or criterion equa- 
tions 1(A)-1(C) will be referred to as 
the “‘end sigma criterion.”’ 

This measure unites in a single concept 
the transient and steady states of the 
system operation as well as the largely 
neglected intermediate states. None 
of these operating states is discriminated 
for or against and because of p(t) each 
gets the weight due it because of the 
relative frequency of its presence at 
such times as the output is utilized. 

In addition to uniting the various 
states of system performance, equation 1 
also unites most of the various criteria 
used to evaluate these states separately. 
As an example, if F(e, t, 1... v,)=e?, 
then referring to equation 1(B) 


s= bi e(t)p(t) dt (2) 


becomes the true mean-square value of 
the error. This end sigma, or prob- 
abilistic error, unifies in one concept 
the mean square criterion used _pre- 
viously for stationary!? and nonsta- 
tionary? random input environments 
and the integral square criterion used 
for nonstationary random? and deter- 
ministic? environments. These become 
special cases of equation 2. Final value 


s= f F(e)g(t)dt 
F 0 


systems also are special cases offequat 
1 by the use of a probability functior 
given in equation 11. However, equatio 
2 is not restricted to these special case 
This is general in that it truly averag 
the square error whatever may be ¢f 
operating state of the system. 

Performance criteria similar in f 
to equation 1, but basically differen 
meaning, have been proposed.7~!# 
basic difference is that instead of 
probability distribution function, { 
of equation 1, an arbitrary weigh 
function is proposed either in a @o 
crete’!4 or a general form.1%/14 ¢ 
sequently these criteria define arbitre 
weighted penalty factors without ¢ 
physical meaning. Nevertheless, e€ 
tions 1(A)—1(C) include as special ¢ 
these arbitrarily weighted penalty 
teria. Assume for instance 


F(e(t),t,21.. .vr) = F(e) WE) 


where W(t) is a function weighting 
penalizing, the errors according to 
time when they occur. W(t)=# 
W(t) =# have been proposed and iny 
gated.’-14_ Then, for instance, eq: te 
1(C) reduces to 


where 
g(t) = W(t)p(t) 


now becomes a modified probabil 
function which, in addition to the relat 
frequency of occurrence of vafi 
utilization times, also includes th 
signers opinion regarding the re 
importance of these utilization t 
Equations 1(C) and 4, being identie 
form, can be handled by ide 
techniques which are introduce 
subsequent parts of this paper. 
significant that these techniques me 
possible to obtain mathematical solt 
for the optimum systems undet 
earlier performance criteria, some 
which has been done heretofore onl 
cut and try. — 

Using equation 1, it is not nece 
to define a number of separate 
formance criteria. Transient and s 
conditions, deterministic or random 
environments, linear and nonlineai 
tems are all covered by the same crit 


Special Cases of Penalty Functi 


If the magnitude of the error al 


then 
F(e)=|e| 


the penalty function, The widely 
1 penalty function 
e (7) 


ycourse puts an increasingly high 
hity on increasing errors. This 
Baty does not fit such cases as 
pmatic machine tools or missiles 
jre all errors up to a certain tolerance 
jalmost equally acceptable and all 
ts beyond the tolerance are equally 
ceptable. For such control systems 


e 


Dh 2 (8) 


i k a constant, would be quite suitable 
e this penalty function (Fig. 1) 
alizes errors progressively up to a 
ain value and then puts an essentially 
tant penalty on errors past this 
e. 

€n more suitable would be the 
wing: 


et 
aa! Ple|tde=pit) (9) 
SH eke 

h gives the probability of error, e, 
g outside the tolerance value, e; at 
#. In equation 9, p(elt) is the 
litional probability distribution in 
nd ¢. For the penalty function of 
ition 9, equation 1 yields 


Pr(t)p(t)dt=1— 
0 


f z if * plelt)p(0)dedt 
2 —erJ0 

a [esaea * (10) 
aw —e} 0 


h is simply the over-all probability 
ne error being outside the tolerance, 
st appropriate criterion, which in the 
of a missile could readily be extended 
e over-all hit probability. 

should be emphasized that the 
ific F(e) function is determined by 
mature of the problem and is to be 
plished by the same engineering 
ment which establishes countless 
- factors of the design. The penalty 
ion F(e)=e? is the most commonly 
and it will be used exclusively in 
emainder of the paper. 


ial Cases of Probability 
istribution Functions 


- definition, p(t) is the probability 
ibution of all times when the out- 
is being utilized, based on a Jarge 
ser of repeated applications of 
is of one given type. The starting 
. of the time scale ¢=0 for deter- 
ig P(t) is the instant of activation 
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eer 


aft 


~I/k 8 1/k 
e 


Fig. 1. Two penalty functions 


of the control system. For a missile, 
‘=0 can signify the instant of firing or 
the instant of transition to the homing 
phase. For some other system, /=0 
may be the instant of turning on the 
power in the control system or of the 
arrival of a step input, etc. 

The probability distribution function, 
p(t), can be readily established in some 
cases. In others, exact determination 
would require an extensive operations 
study, but a reasonable estimate can 
almost always be found quite readily. 
As a matter of fact, by a slight shift of 
the sense of the definition, p(t) could be 
represented as expressing the distribution 
to which the spread of the possible 
utilization times can be narrowed down on 
the basis of the limited information 
available to the designer. In this sense, 
P(t) would become, at least partly, 
representative of the lack of definition 
in the design parameters rather than a 
probability distribution of the actual 
operating times. - 

As an example, consider a control 
system which receives step inputs at 


t=0 followed by one instantaneous 
measurement at f= T». 

Then 
p(t) =4(t— Ty) (11) 


where 6(¢) denotes the Dirac delta func- 
tion. If in the same system the measure- 
ment extends from 7) to J; (as it would 
for instance in some wind-tunnel appli- 
cations) 


P(t) = [u(t— To) —u(t—T1)] (12) 


Tes d'o 
where u(t) denotes the unit-step function. 
In writing equations 11 and 12 the 
assumption was that the times of utiliza- 
tion are the same for every repeated 
application of the step input. Should 
these be variable in some way (in the 
case of a missile the instant of explosion, 
To in equation 11, would vary from 
launching to launching depending on the 
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distance of the target at the time of 
launching) then equations 11 and 12 
modify respectively into 


K 


uf ——— ee 
p(t) = lim =) 8(t— Tox) =6 (¢— 1) 
iy JEG 
7 
(13) 


K 
p(n) = Ii >» el 
Bree Tyx4—Tox 

= 


[u(t— Tox) —u(t—Tx.)] 


[u(t— Tox) —u(t—Tox)] (14) 


ts 1 

Tip Top 
For a missile where maximum and 
minimum ranges (7; and 7) respectively) 
are clearly set, the results of equation 
13 would come close to the form of 
equation 12. 

For systems which operate essentially 

in the steady state, one possible form can 
be produced from equation 12: 


u(t)—u(t—T))] (15) 


1 
j= lim — 
p(t) ere 


or from a Rayleigh probability distribu- 
tion: 


t t 
SSibbiny = ee t>0 
fi tin, 5 ex ( a8) 
=0 t<0 


(16) 


With either equation 15 or 16, equation 
2 goes over into the customary definition 
of the rms error as used in the study of 
systems subjected to stationary random 
processes. This is especially obvious in 
connection with equation 15. 


Finding the Probabilistic Error 
for Linear Systems in s Domain 


The end-sigma mean-square error is 
defined by equation 2. Its computation 
for linear-time-invariant lumped-con- 
stant systems can be carried out either 
in the s domain of the Laplace trans- 
forms or in the time domain. The s- 
domain approach will be discussed first: 

The Laplace transform of the error 
can be expressed as 

M(s) 


E(s)=I(s)—K(s)R(s)= =~ 


Ns) CD) 


where M(s) and N(s) are entire functions 
with M(s) having a higher order singu- 
larity than M(s). I(s) and R(s) are, 
respectively, the Laplace transforms of 
the desired output and the actual input, 
and K(s) is the closed-loop transfer 
function of the system. Now, as shown 
in the Appendix, equation 2 for a deter- 
ministic input [e?(t)=e2(t)] can be con- 
verted by complex convolution into 
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s=lim £2~1{ P(—s)[E(s)*E(s)] j 


1=0 


(18) 


where £2! denotes the inverse 2-sided 
Laplace transform.% P(s) is the 1-sided 
Laplace transform of /(#), and the 
asterisk denotes complex convolution. 
Since p(é) is zero for negative time 
(t=0 being the instant of activation, 
so that any utilization of the output can 
come only in positive time), 


P(—s)=L2[p( —#)] 
A(s) 


rs 3 :: st 
BCS plije*‘dt (19) 


where A(s) and B(s) are entire functions. 
Further, 


1 C+J © 
Bs Ble) = 5 f E(s—w)E(w)dw 
C—jo 


DOOM —0)MO) _ Hs) 
2 UNNI Lis) ae 


where J(s) and L(s) are entire functions, 
the ; are the zeros of M(s), all of which 
are assumed to be single, and 

CS) = = Ns) (21) 
By using inverse transformation through 
contour integration, equation 18 goes 
over into one of the two alternate forms of 
equation 22 and equation 24 


Age) J(gx) 
B'(qx) L( qx) 


(22) 


where the g, are the poles of P(—s), 
all of which are assumed to be single, 
and all are located in the right half plane. 
Also 


Bs) =— £ Bs) (23) 
Alternatively, 
A(pi) Sb) (24) 


B(p1) L'(p3) 


where the p; are the zeros of L(s), all 
of which are assumed single, and located 
in the left half plane. Also 
L"s)= © L(s) (25) 
ds 

Two equations (equations 22 and 24) 
result from the 2-sided inverse trans- 
formation since the function is continuous 
at #=0 and, consequently, its limit 
can be evaluated either from the trans- 
form for positive time or from the trans- 
form for negative time, in other words, by 
closing the contour in the s plane on 
either the right or left. 

Equations 22 and 24 represent the 
general s-plane formulas for evaluating 
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the end-sigma error for time-invariant 
lumped-constant linear systems and de- 
terministic inputs. The evaluation neces- 
sitates the finding of the poles of the 
error E(s) and of the transform of the 
probability distribution function P(—s). 
These are then substituted into the 
rational function-type expressions of 
equations 22 or 24. The evaluation of 
the poles of E(s) usually can be ac- 
complished by root-locus techniques. 
The substitution into equations 22 or 24 
can be aided materially by conventional 
spirule operations. Inherent in these 
derivations is the assumption that E(s) 
has no pole at the origin. 


For example, assume an input 
r(t) =a(t) = 14 2¢ t>0 


and probability function 
pi)= ete") 


while the system closed-loop transfer 


function is 


Sarl 
ei Leeeren 
Then 
RG) —=T6)'— Pd 
2 
E(s)=1(s)—K(s)R(s)= 3 2 arty 


so that p:=—0. 5+jV0.75 (5 pe=—0.6— 
jV0.75 


Also 


10 
(s—2)(s—5) 
so that m= +2, q@=+5. 


P(-s)= 


From equation 20 
(s —2p2)(s—pit2)(pi+2)— 


I(s) _ (s—2p1)(s —p2+2)(p2+2) 
L(s) (s—2p1)(s—2p2)(s —p1—p2)(p1— pr) 


and the end-sigma error from equation 
22 is s=1.37. 

This highly simplified example should 
be interpreted as illustrative of the 
computation process and not as illustra- 
tive of the power of the method attainable 
when working equations are developed 
from general solutions given in this 
paper. 


Finding the Value of the 
Probabilistic Error in the t 
Domain 


If the end-sigma error defined by 
equation 2 in the ¢ domain is to be 


the input, and k(#)=0 for t<0O whict 


_ evaluation of the integrals in equa 


duration, (p(f)=0 if t>T), by st 


evaluated in the time domain, first j 
will be assumed that the system is tit 
invariant linear with known impuls 
response or weighting function, & 
Consequently the error can be defiz 
by time-domain convolution, or suf 
position, as 


dai [ R(ti)r(t—h)dty 
0 ‘ 
where i(t) is the desired output, 


necessary for a realizable k(t). 
Let the following function be defined 


eae f ” lta tel tee tet vate 
: 


Then 
s=0%(0,0) 


Substituting equation 26 into equat 
27, and rearranging by an exchange 
the sequence of integration and con: 
ering equation 28, there results 


s=Wii(0,0)—2 f dh k(t) ¥ir(0,—h) + 
0 


{ and) { di,k(t; \Wrr(tsst3— by) 
0 0 


where the following definitions . 
made: 


vottrd= f dtp te+t)r(ta)r(te-+ 7) 
vadtr)= | auplictDnOKeD 
0 


Wirlt,7) = if dizp(ta+ t)i(te)r(to+ rH 


where the straight line denotes ense! 
averaging for the conditions of eque 
1(B). 

Equation 29, then, is the time-dot 
equivalent of equations 22 or 24, 


29 yields the end-sigma error value f 
given system. The results are ree 
generalized to cases of finite m 
transfer functions, (k(t)=0 if t>T 
probability density functions of f 


C} 


changing the upper limit of the inte 
to T in equations 1 and 27-32. 

Some additional discussion, how 
is due equations 30-32. These equé 
define a generalized form of auto 
cross-correlation functions which 
in one concept the auto- and 
correlation functions customarily us¢ 
stationary processes! and the 
and cross-translation functions prof 
for transient conditions.? Ay 
probability function of ‘equation 


nediately reduces equations 30 through 
o the definitions of auto- and cross- 
elation functions. The same prob- 
ty functions would result in a 
ified but equivalent definition of the 
islation functions.? An infinitude of 
i special cases result from the various 
Hicular probability functions. 

aking a time varying Laplace trans- 
p* of equations 30 through 32, a 
alized power-spectral density cor- 
jonding to the general correlation 
tions would result. 

eferring to equation 30, 


| 


5S) = f Vrr(t, r)e**dr 


= IT  p(ti-+t)®(n,s)dt (33) 
re 
“S { __ ticker (34) 
he apparent power-spectral density 
the ensemble at time 4. If the 


mble is stationary and the ergodic 
ry applies, 


(35) 


the generalized power spectral 
ity, 


as) f ” plhte)dn=4(s) (36) 


ices to the conventional power- 
ral density of the stationary process, 
of course should. Similar definitions 
easily obtained from equation 32 
a generalized cross-power spectral 
sity. 
g the Optimum Invariant 
near System Under the End- 
gma Criterion 


quations 22, 24, and 29 are useful in 
lating the end-sigma error for a 
m of known construction and charac- 
tics by working respectively in the 
tdomains. If a choice exists between 
ral potential systems, finding and 
paring s for each of these immediately 
pinpoint that particular system 
+h gives optimum error performance, 
is, minimum true rms error. This 
ld be an optimum from among a 
ted choice, established by cut-and- 
techniques. Now if the choice of 
sms is made unlimited to the point 
meluding all linear systems with 
tant coefficients, then a direct solu- 
for the optimum system becomes 
ible, by finding the optimum A(t) 
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weighting function or its transform, the 
K(s) transfer function. Since either of 
these unknown functions is included in 
the defining integral of s, the quantity 
to be minimized (equations 18 and 29 
respectively), the solution can be ob- 
tained by calculus of variations. This 
will be illustrated here for the time 
domain, equation 29. The minimum of 
s is defined by equating its first variation, 
és, to zero, that is from equation 29 


ds=2 { dtsk(t) X 
0 


ae ani Weft —1)—vel0,—0 | =0 
0 


(37) 


which is only possible for arbitrary 
dR(t) if 


ib dihk(ti)Wrr(t.té—t1) —Wir(0,-t)=0 t>0 
0 
(38) 


It can be shown that equation 38 
actually defines a minimum. 

As in the case of equation 29, equation 
38 is also easily adapted to finite memory- 
transfer functions and finite duration 
probability distribution. 

Equation 38 is a linear integral equation 
of the first kind. ¥,,;(¢,t—#) is its kernel 
and as may readily be shown from equa- 
tion 30 this kernel is symmetric. The 
solution of this integral equation for 
k(t) is the impulsive response of the 
optimum system, the system yielding the 
minimum end-sigma error. The latter 
then should be found from equation 29. 

The literature, both mathematical! 
and engineering,* for solving integral 
equations is extensive. No attempt 
toward completeness is made in the 
references, 


Conclusions 


As a measure of performance, the end 
sigma or probabilistic error eliminates 
the indirectness so characteristic of usual 
design techniques. The end-sigma error 
summarizes the entire system perform- 
ance, rather than explaining it in full 
detail, thus making it a desirable starting 
point for analysis or design. Such 
techniques as s-plane transient analysis 
can then be used to spotlight special 
limited aspects of the performance, 

In addition to introducing this measure 
of performance it has been demonstrated 
that it can be evaluated in either the s 
or ¢ domain for time-invariant linear 
systems. An integral equation of a 
common type has been introduced to 
define the impulsive response of the 


Zaborszky, Diesel—Measure of Control-System Performance 


optimum-invariant linear system under 
the end-sigma criterion. 


Appendix 
Equation 2 can be rephrased as 


s=lim | er-per nar (39) 
¢=0 0 


Or introducing 2-sided Laplace transforma- 
tion 17 (o[ ] and inverse 2-sided Laplace 
transformation £.—1[ ] 


s= him Gs," ol f Mobeni 
t=0 0 
= lim cn} f e *“dixX 
t=0 —o 
f Heber nar] 
0 


= lim el f e(7)e—**drX 
t=0 0 


le r= neto4te—n | (40) 


—co 


Further, from complex convolution theory 
ealern) [ e(r)e*"dr=E(s)*E(s) (41) 
0 


Equations 40, 41, and 19 can be combined 
into equation 18. 
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Modern Systems of Traffic Control as 
Applied to the Seaboard Air Line 


Railroad Company 


J. R. DEPRIEST 


NONMEMBER AIEE 


N December 6, 1941, the Seaboard 

Air Line Railroad placed in service 
the first segment of a major installation 
of a system of traffic control on its 
line. The segment was only 33 miles in 
length but it made history. It was a 
portion of a 65-mile segment of railroad 
to be electrically controlled from a 
remote location 157 miles from the most 
distant point. At that time this line 
became the longest control line for this 
type of service in the world. The system 
proved to be very efficient and econom- 
ical, and expansion of the system was 
made as rapidly as material and labor 
could be obtained. As of today, trains 
operate the entire distance between 
Richmond and Miami-Tampa by signal 
indication controlled by modern systems 
of traffic control. 

Prior to December 6, 1941, the Sea- 
board line consisted of a single-track 
railroad with more than 4,000 miles of 
track extending from Richmond to 
Miami and Tampa. A main line also 
extended from Hamlet, N. C., toward 
Atlanta and Birmingham together with 
other extensions operating in the states 
‘of Virginia, North Carolina, South Caro- 
lina, Georgia, Alabama, and Florida. 
At that time the Seaboard had approx- 
imately 550 miles of automatic signals 
mainly extending between Richmond, 
Va., and Monroe, N. C.; between 
Savannah, Ga. and Baldwin, Fla.; 
and between Starke and Coleman, Fla.; 
see Fig. 1. With the systems then in 
service it was necessary for telegraph 
operators or telephone operators to be 
placed at strategic points along the rail- 
road in order to receive train orders 
issued by a train dispatcher. The 
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operator then physically delivered these 
orders to train crews in order that op- 
posing trains could meet at strategically 
placed sidings at the proper time. Since 
it takes considerable time to issue and 
check train orders, it can be seen that 
the system becomes difficult and ex- 
pensive to operate as the number of 
trains increases. This is so because the 
total number of meets in a fixed distance 
increases directly in proportion to the 
square of the number of trains. Perhaps, 
this point should be further discussed 
so that the reader may become familiar 
with the problems encountered in de- 
signing a traffic-control-system for op- 
erating a single-track railroad. 

It is evident that for a double-track 
railroad it is possible to operate trains of 
equal speeds at eqttal time intervals 
without interference from either fol- 
lowing or opposing trains, since trains 
moving in opposite direction operate on 
adjacent tracks. Fig. 2(A) is a time- 
distance chart showing the position of 
southward (1S) and northward (1) 
trains. Train 1S leaves station A 
proceeding toward station C at a constant 
tate of speed. The progress of this 
train-is represented by a diagonal line 
with a negative slope. At the same time 
train 1N, represented by a diagonal 
line with a positive slope, operating in 
the opposite direction, leaves station C 
proceeding toward station A at the same 
speed. Since each train is operating at 
the same constant speed they will arrive 
at their destinations C and A, respec- 
tively, after a time interval ¢, having 
traversed the distance D between sta- 
tions. It will be noted from Fig, 2(A) 
that the two diagonal lines, representing 
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the progress of train 1.5 and 1N, interse 
or cross one another half way bety 
stations A and C, after one half of 
time interval ¢ has elapsed. Since f 
1S leaving station A will meet only 
train between stations A and C (and 
train at C if the railroad continues f 
greater length) the railroad could 
single track with the exception ¢ 
small distance near the center (siding 
and at stations A and C if means ¥ 
provided whereby one of the trains e¢ 
enter, pass the opposing train, and 
leave the siding. It should be poi 
out that if train 1S continued bey 
station C another train would be me 
station C if opposing trains were 
patched at equal intervals. ; 

Similarly Fig. 2(B) illustrates 
conditions encountered if the numbe 
trains operated in time interval ¢ y 
doubled. It can be seen from Fig. 2 
that the total number of meets occur 
within a time interval ¢ would iner 
from two to eight. However, 
sidings would be required as comp 
with two sidings previously. 
trains in each direction were oper 
within time interval ¢, then by refer. 
to Fig. 2(C) it can be seen that a tot 
18 meets would occur within time inte 
tand 6 sidings would be needed. 

It is interesting to note that where 
train in each direction is operati 
time interval 7, siding B is used one 
as shown by the heavy dot loc 
vertically under siding B. Where 
trains in each direction are ope 
per time interval ¢, siding B is 
twice. And where three trains 
operating in -each direction per 
interval ¢, three trains use siding B. 

If ¢ is made equal to a time interv: 
one day and Table I is prepared fi 
the foregoing information, a fort 
can be written from the first three ¢ 
tions to determine the total numbs 
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S occurring within the area the 
2 operates in one day or for any 
interval selected. Similarly, for- 
aS can be derived showing the 
ts per train per day which equals the 
ber of sidings needed within the 
These formulas are as follows: 


where 


t=time interval a train operates in traveling 
distance D 

T=number of trains operated per direction 
per time interval ¢ 

M=total number of meets occurring within 
time interval ¢ or in distance D 

D=distance train travels in time interval ¢ 

S=number of sidings needed within dis- 


(Tt)? tance D or time interval t 
prt For a practical example, assume ¢ 
Tt is equal to one day and one train is 
ION STATION STATION STATION 
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Fig. 2. Time-distance chart 
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operated each half-hour, then 48 trains 
will be operated in each direction per 
day. From Table I, under S, 96 sidings 
will be needed. If the trains average 40 
mph, a train will travel 960 miles in one 
day, which is distance D. Therefore, if 
the 960 miles is divided by the 96 sidings 
it will be found that sidings must be 
spaced 10 miles apart. It should be 
pointed out that with sidings spaced 10 
miles apart and trains averaging 40 
mph, if one of the trains was delayed and 
could not make the scheduled progress, 
it would be necessary to require the 
delayed train to enter the siding one 
station before its scheduled point to 
avoid delay to the opposing train. 
This train would then be further delayed 
for a time interval not exceeding the 
sum of its running time to the next 
station and the running time of the 
opposing train from the next station to 
the siding where the meet will occur. 
Where the speeds of each train are 
equal, then the maximum additional 
delay that could occur would be two 
times the running time between sidings, 
or in this case 18 minutes. To reduce 
this type of delay more sidings are 
usually provided than required to handle 
the maximum number of trains as shown 
by Table I. 


STATION STATION 


L= TIME 


fD=DISTANCE A TRAIN TRAVELS PER TIME INTERVAL td 
() 
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A traffic-control system remotely con- 
trolled must provide electric power to 
operate and control a railroad switch 
located at the entrance or exit of a 
siding, and extend information to en- 
ginemen of trains by means of wayside 
signals appropriately located. It is 
necessary for the dispatcher to control 


the following functions for each siding 
switch: 


1. Position switch normal. 
2. Position switch reverse. 


3. Position 
stop or clear. 


northward signals to either 


4. Position 
stop or clear. 


southward signals to either 


5. Position automatic reclearing feature 
for following trains. 


6. Operate maintainer’s call light. 


The system must indicate certain 
information from the field location, 
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Fig. 3. Track diagram with 
associated indication lights 
and control switches 


where the power-operated switch is 
located, to the dispatcher located at 
control machine at a remote location. 
The following information is required 
from each power-operated siding switch: 


1. Occupancy or nonoccupancy by train 
located on track north of switch. 


2. Occupancy or nonoccupancy by train 
located on switch. 


83. Occupancy or nonoccupancy of train 
located south of switch either on main line or 
siding. (Usually a typical procedure is to 
assign siding occupancy indications to 
apparatus located at north switch of siding 
and occupancy of track between switches to 
apparatus located at south switch. ) 


4. Indication to show when switch is set 
for main-track traffic or in transit position. 


5. Indication to show when switch is set 
for side-track traffic or in transit posi- 
tion. i 


6. Indication of northward signal aspect 
(clear or stop). 
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SECTION E————- 


126 MILES 
INDIANTOWN DYER 


7. Indication of southward signal as 
(clear or stop). _ 


8. Power off indication (on or off). 


Fig. 3 illustrates a track diagram 
associated indication lights and c 
switches. a 

From the foregoing it will be 
that 6 bits of information are ne 
for control purposes and 8 bits — 
formation are needed for indi 
purposes at each siding switch 
in traffic-control territory. The 
Switch 506-A system, describe 
Mr. Baughman! in another pape 
arranged to control and indicate t 
7 bits of information for 35 station 
one pair of wires. From the for 
requirements it can be seen that 
the exception of indications, this 
will be adequate. Since power is 
small percentage of time, it is pr 
to indicate both a normal and 
switch condition simultaneously, 
since this condition is physically im 
sible a means is established to ¢ 
that the a-c power is off. Since 
system uses a share-time basis for ¢ 
trols and indications, it is usually 
practice not to control more th 
switches from. one line, or 15 sic 
Experience has proved that this nu 
of sidings can be controlled 1 
unduly delaying indication or 
codes for all*stations on a system 
by a busy railroad. If more 
sidings are needed, the territon 
divided in two or more sections 
a system is applied to each see 
When this is done, a pair of wire 
tending from the most distant si 
switch to the control point is nec 
for each section of line to which a s 
is applied. If the d-c coded in 
at the end of each section are 


MIAMLFLA. 


| TYPE "co" | 
REPEATER 


| 
| 
! 
| 
| 


Fig. 5. Traffic control machines located in Jacksonville 


llate a voice-frequency telegraph- 
carrier, only one pair of wires 
be needed throughout the entire 
vad for this purpose. For control 
, similar carrier apparatus is used 
the modulation takes place at the 
ol machine and the carrier is 
dulated at an end of the remote 
on in the field. Fig. 4 illustrates the 
ol arrangement used by the Sea- 
1 Air Line Railroad between Jack- 
He and Miami, Fla. 
will be noted in the foregoing ar- 
ment that a Western Electric C 
which provides three trunk tele- 
© circuits was applied to the single 
bf wires of the code line and cable 
ding between Jacksonville and Wild- 
, Wildwood and Sherman, and 
an to Miami. One channel of 
atrier was used to apply the voice- 
ency signal-control carriers. The 
line or system called AA extends 
een Jacksonville and Baldwin. The 
narked A extends between Baldwin 
Hawthorne with carrier controls 
ding this line into Jacksonville 
its north end at Baldwin. The 
line designated B extends between 
horne and Wildwood. Since the 
‘tier terminal is located at Wild- 
_ the south end of this line is ex- 
d to Jacksonville by carrier. The 
esignated C extends between Wild- 
and Winter Haven. This line is 
led from the north end by carrier 
Jacksonville. Line D_ extends 
en Winter Haven and Sherman and 
eatrier C terminal is located at 
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Sherman, it is extended from the south 
end to Jacksonville by carrier. The 
next line designated E extends between 
Sherman and Yamato and the north 
end is extended by carrier from Sherman 
to Jacksonville. The last line designated 
F extends between Yamato and Miami 
and since the C carrier terminal is 
located at Miami the F is extended from 
its south end to Jacksonville. The 
distance from Jacksonville to Miami is 
approximately 405 miles. The controls 
for this entire territory are placed on 
three machines, one controlling the terri- 
tory between Jacksonville and Baldwin, 
another between Baldwin and Wild- 
wood, and the other controlling the 
territory between Wildwood and Miami. 
Fig. 5 shows the traffic control machines 
located in Jacksonville. 

Up to this point the communications 
involved have been discussed, but not 


the actual wayside signal systems avail- 


able for traffic control. On the Seaboard 
three distinct types of systems have been 
installed. The first system is referred to 
as the full-centralized traffic control 


_system which permits trains to follow one 


another without waiting until the pre- 
ceding train has reached the next siding. 
This system is equipped with power- 
operated switches and, where required, 
the sidings are usually more closely 
spaced than in the other systems de- 
scribed later. The next system has been 
designated as the modified-centralized 
traffic control system. This system 
provides for power-operated switches, 
but will not permit trains to follow until 


the preceding train has reached and 
passed the entrance switch to the next 
siding. The system eliminates the many 
wayside signals needed when following 
movements are permitted as in the full- 
centralized traffic control. The third 
type of centralized traffic control is 
designated by the Seaboard as remote- 
controlled automatic signals system. 
This system basically is arranged for 
train operation by signal indication in 
the same manner as the modified- 
centralized traffic-control system, except 
that following trains cannot follow 
another train until the first train has 
reached and passed the north end of 
the siding ahead. It can be seen that 
this system has eliminated much appara- 
tus at one end of each siding. The 
siding switches are usually hand operated. 
However, the system is so designed that 
power-operated switches can be added 
when needed at a minimum of expense 
and the system can be expanded into 
either modified-centralized traffic con- 
trol or full-centralized traffic control by 
adding the necessary apparatus for such 
system. In addition to the traffic control 
systems the Seaboard has short stretches 
of double track where the conventional 
automatic signals are provided and where 
trains operate by signal indication with 
current of traffic. Fig. 6 shows the 
various types of centralized traffic con- 
trol used on the Seaboard as of October 
1958. In addition the map shows the 
absolute-permissive block signal system 
between Hamlet and Atlanta. This 
system is a conventional automatic- 
permissive block signal system and is 
used with train orders for the movement 
of trains. It will be noted from this map 
that the Seaboard operates by signal 
indication for the entire distance of 
more than 1,040 miles between Rich- 
mond, Miami, and Tampa. 


Three control machines are located 
at Raleigh, N. C., which control the 
entire territory between Richmond, Va., 
and Hamlet, N. C. Three control 
machines are located at Savannah, Ga., 
which control the territory between 
Hamlet via Columbia to Savannah; 
Hamlet via Charleston to Savannah; 
Savannah-Gross—Jacksonville and Gross- 
Baldwin, respectively. At Jacksonville 
there are three traffic control machines, 
which control the entire territory be- 
tween Jacksonville and Miami. One 
machine controls the Jacksonville-Bald- 
win section, one the section between 
Baldwin and Wildwood and the other the 
section between Wildwood and Miami. 
The section between Coleman and 
Tampa, Florida, is controlled from a 
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machine located in Tampa. Incidentally, 
at the present time the control line 
extending between Jacksonville and 
Miami is the longest centralized traffic- 
control line in the world to extend in one 
direction from the control point. Also 
the system between Richmond and 
Miami is the longest continuous system 
in the world on which trains are operated 
by signal indication. 

After the war, in 1946, when materials 
were again available the Seaboard under- 
took the installation of traffic control 
on many of its main lines where signals 
were not in service. The territories to 
be controlled were long and a reliable 
code line was needed. Methods were 
developed to economically install an 
underground cable consisting of two no. 
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Fig. 6. Various 
types of centralized 
traffic control used 
on Seaboard as of 
Oct. 1958 
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10 American Wire Gauge solid copper 
wires, each circumscribed by 3/64-inch 
polyethylene insulation. The two poly- 
ethylene insulated conductors were 
twisted with a 4/64-inch Flamenol jacket 
covering the twisted pair. Details of 
this installation were given in another 
paper.? Part of this cable has been in 
service for 11 years, and its performance 
has been satisfactory with the exception 
of that cable installed in certain por- 
tions of Florida where the ground was 
infested with pocket gophers. In this 
territory the pocket gophers, which nor- 
mally feed on roots, attacked and 
ate the cable damaging the insulation 
and eventually causing grounds. The 
resultant electrolysis eroded the con- 
ductors. Experience has proved that 


‘is clear for either following or opp 


the same cable with a 10-mil bre 
tape wrapped spirally around the. cak 


cable from pocket gophers. 
has been found that pocket gop 
gnawed the armored cable but ~ 
unable to penetrate the armor. O;) 
the armored cable installed onl 
failure has occurred and that as a resu 
a nail having penetrated the armor 
insulation, evidently while on the 
before being installed. It is interes 
to note that no lightning failures on 
armored cable have occurred but 
the nonarmored cable an occas} 
failure does occur due to lightning, 
though a test sample of the cable ¢ 
6 years of service was examined 
found to, have a breakdown voltag 
excess of 40,000 volts to ground. 


It should be pointed out that 
systems of centralized traffic control 
wayside signal apparatus used is of 
quality. The electric circuits are 
signed on a “‘fail safe’’ basis and, — 
though the communication facilities 
trolling such wayside signs may fa 
operate properly, the wayside signal 
arranged to check that all power-ope: 
switches are mechanically locked 
hand-operated switches are in the pi 
position, rails are intact, the track a 


trains, and that the signals gove 
opposing movements have been at 
for the prescribed predetermined © 
In the traffic-control systems inst 
recently on the Seaboard, in the € 
of communication failure the — 
automatically will clear signal in adi 
if the route ahead is clear. This fez 
is very desirable and eliminates a 
number of delays if a communie; 
failure should occur and it meets 
company’s first and last requisite w 
is safety. 
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sification yards for the handling of freight 
cars. These installations include meth- 
ods for switching each car to its proper 
destination track in response to previously 
stored information. This storage may be 
on teletype tape or other suitable memory 
system. Radar devices are used for speed 
measurement in servo loops for the con- 
trol of the speeds of the cars. These in- 
stallations are among the nearest ap- 
proaches to full automation in service 
today. 


Description and Discussion 


The preceding indicated some of the 
electrical developments that have been 
adopted and developed by the railroads. 
The paper shall now discuss centralized 
traffic control (CTC). This is actually a 
system of train operation where the train 
crew receives no written orders from the 
dispatcher. The train is operated in ac- 
cordance with the signal indications, 
Train orders may be considered as issued 
by the dispatcher because he is given 
remote control of the signals and switches 
in his territory. A basic concept of this 
system is that it is divded into two main 
parts. The first is a code communication 
system so the dispatcher may convey his 
desires to the signals and switches. Also 
this communication system transmits in- 
formation to the dispatcher giving the 
positions of signals and switches as well 
as the location of all trains. The second 


part consists of a complete wayside signal 
system with the main signals normally in 
the “stop” position. The over-all safety 
resides to a major extent in this part of the 
system. There are two conditions that 
must be satisfied before a signal can clear 
to authorize the passage of a train. The 
dispatcher must send a code to the par- 
ticular signal, in effect requesting it to 
clear. This request is received and de- 
coded automatically and thus the first 
condition is satisfied. The second condi- 
tion involves the automatic signal equip- 
ment which functions continuously in re- 
sponse to traffic conditions. Before the 
particular signal can clear, the automatic 
signaling equipment must check that it 
is safe for the signal to clear. In a simi- 
lar manner the dispatcher has control of 
the track switches in his territory. But 
before any track switch changes its posi- 
tion, in response to a code from the dis- 
patcher, the automatic signaling equip- 
ment must check that it is safe for the 
switch to move. 


Fig. 1 shows a machine at Savannah, 
Ga., for the control of 250 miles of single- 
track railroad, A diagram on the machine 
shows the track layout over which the 
dispatcher has control. Small lights on 
the diagram show the positions of all 
trains and these positions are recorded on 
a chart usually called an automatic train 
graph. Also the positions of all switches 
and signals are indicated by lights on the 
panel. 


The dispatcher sends codes to the 
various locations in his control of the 
traffic. This is accomplished merely by 
the operation of small levers and push 
buttons. Control codes are generated 
automatically in response to the positions 
of the levers and push buttons. Indica- 


Fig. 1. A 25-foot control machine in service at Savannah, Ga. 
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Fig. 2. A multiple CTC line circuit with low-pass filter to permit a physical voice circuit on the 
: same pair of line wires 


tion codes from the way stations to the 
control office are generated in response to 
the positions of the functions in the field. 

Two line wires are used for the trans- 
mission of these codes. The use of a 
single wire with ground return received 
consideration but was discarded promptly 
because of the requirement of co-ordina- 
tion with communications lines. Another 
reason against the single wire with ground 
return is the interference that would be 
experienced during magnetic storms. 

The usual installation on a single-track 
railroad has the field functions arranged 
in groups of about 6 to 10 functions per 
station. The stations are spaced along 
the line at intervals in the order of 5 
miles. In usual communications lan- 
guage this might be called a party line. 
In this application there are frequently 15 
to 25 parties on a line. 

Fig. 2 shows a simplified line circuit. 
In the event of simultaneous starts at 
office and field, the office is given prefer- 
ence by making the first pulse of each 
control code long as compared with the 
first pulse of each indication code. The 
method for establishing priority among 
the field stations is dependent upon the 
specific assignment of station codes and 
the geographical sequence of the stations 
with respect to the office. 

The codes consist of pulses of long and 
short intervals. The first portion of a 
code is assigned to station selection and 
the last portion to the functions of that 
station. 

The main novelty of Fig. 2 with re- 
spect to previous line circuits is that the 
way stations are connected to the line in 
multiple. The two line wires are con- 
tinuous throughout the installation whe- 
ther the line is idle or handling codes. 
Way stations transmit indication codes 
by alternately shunting and removing the 
shunt from the line. This action produces 
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changes in current from the line battery 
at the office. These changes in current 
through the indication transformer in- 
duce voltages in its relay winding for 
operation of the office line relay. 

Becatise of the over-all reliability re- 
quired of this system special attention is 
given to the strength and ruggedness of 
the line wires. This being true it is logical 
to consider the use of these same line 
wires for other functions. This principle 
develops in a cumulative manner. 
stronger and more rugged the line wires, 
the greater is the desirability of placing 
additional functions on that same pair of 
wires. Also the greater the number of 
functions being handled over the one pair 


of wires, the greater is the attention that — 


is warranted in providing reliability and 
good maintenance. 

A major step in permitting additional 
functions to be handled on the same pair 
of line wires is the installation of low- 
pass filters between each way station 
equipment and the line. The low-pass 
filter for this service has low resistance to 


_the flow of direct current and high imped- 
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ance to frequencies in the voice and ultr 
sonic ranges. The filter is designed 
minimize the electrical noise which is f 
to the line because of the coding contagt 
To maintain ruggedness and reliabilit 
of the line circuit, special attentiot 
given to this filter so that it will withst 
the severe surges that develop on line 
this type during storms. This special 
tention to the filter does much to pro 
other apparatus at the station fro 
damage by lightning surges. With t 
low-pass filters in place, as shown in 
2, the pair of wires for the CTC line se 
as an excellent voice communications 
cuit. It is the usual practice to ha 
telephone at each way station so 
maintenance men and train crews 
communicate directly with the dispate 

The system as described up to- 
point was available in 1940, largely 
cause of the intelligent use of funds 
personnel for intensive development 1 
during the depression years of the 1 
ties. During the early days of World 
Il, the traffic congestion on the st 
track railroads was so great that Ir 
studies were made to determine the 
way to relieve the congestion prom 
and consistent with the economic u 
strategic materials. The traffic cot 
system was determined as the most f 
Emergency measures 
much overtime work by the man 
turers and the installation crews pro} 
relief at the most critical locations. 

The system operated satisfactoril 
practice, but soon after there was ¢ 
mand for further expansion beyon¢ 
approximately 75-mile limit of the sy 
as it existed. The experience indic 
that one dispatcher could handle 1 
more territory and that additional a¢ 
tages would be realized with larger 
centration of controls. ; ; 

The solution to permit large ins 
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4. Total road miles of traffic control 
Versus years 


ns, hundreds of miles in length, to be 
rated from one control office gave full 
sideration to the desirability of retain- 
as standard the system and equipment 
t had been giving excellent service. 
er a great deal of study of various ap- 
aches to a solution, the one that 
ved to be most acceptable considered 
system as it then existed as a unit of a 


to end and controlled from one office. 
= hear unit operates in the conventional 
nner, previously described. The sec- 
unit is arranged in exactly the same 
nner as if there were a control point 
ted at the far end of the first unit. 
p carrier frequency transmitted in 
h direction over the first unit permits 
second unit to be controlled over the 
wires of the first. In a similar man- 
a pair of carrier frequencies over the 
Wires of the first and second units 
its the handling of a third unit. 
filters which were provided to per- 
a telephone communication system 
be operated on the one pair of wires 
serve to permit the transmission of 
carrier frequencies. This use of car- 
frequencies removed the previous 
ations as to size of the installations. 
roads may now handle any desired 
of installation from one control point. 
outstanding example is the Savannah 
e on the Seaboard Air Line railroad 
1 which over 600 miles of track are 
rolled. 
9 retain the telephone circuit on the 
e pair of line wires it was necessary to 
t the frequencies of the carriers so 
they would not be audible to the ex- 
of causing interference or annoy- 
to telephone conversations. Tests 
experience indicated 11 kc to be the 
mum frequency that was practical. 
interesting to note that the first car- 


ry 1959 


rier units for this type of service were 
built on a temporary basis to serve until 
more substantial units could be made. 
These units were installed on the Southern 
Pacific in 1942 and are reported to be 
continuing in satisfactory service at the 
present time. To minimize the delay to 
trains that might be caused by the failure 
of electronic components the carrier units 
were installed in duplicate. 

Starting with the first carrier units 
there was an attempt to simplify the 
manufacturing and inventory problems 
by making as large a portion of the unit 
as possible common to all frequencies, 
The frequency determining units were 
made separately as plug-in devices. 

Special attention should be given to the 
way that the use of carrier frequencies ex- 
panded the size and scope of CTC. All 
of the main items retained the same use as 
in the previous system. The addition of 
carrier frequency and coordinating equip- 
ment broke through the previous barriers 
of the size of installation and the distance 
of transmission. This indicates why it 
was possible to meet the sudden demands 
for larger installations. An idea of the 
sudden change in demand is indicated by 
noting the change in slope of the curve in 
Fig. 4 during the years 1942 to 1943. 

An important application of the equip- 
ment and system just described, outside 
the railroad field, was in the extensive re- 
mote control and indication system for 
some of the vital experiments in the Man- 
hattan Project. 

An interesting demonstration was per- 
formed jointly by the Pennsylvania Rail- 
road, the Radio Corporation of America, 
Western Union, and Union Switch & 
Signal in 1946 to show that the control 


office could be separated from the active 
territory by great distances and various 
forms of communications. The transmis- 
sion path between the office and the ac- 
tive territory in the demonstration in- 
cluded carrier frequencies on railway line 
wires, Western Union telegraph lines, and 
several sections of microwave. 

The 90-mile CTC installation which 
was made on the Orinoco Project in Vene- 
zuela, with the co-operation of the Radio 
Corporation of America, in 1953, is novel 
with respect to all other installations be- 
cause it uses microwave principles instead 
of line wires. 

A recent and significant development 
has been made in the control machine to 
facilitate the work of the dispatcher. 
Fig. 5 shows this new concept of machine 
which has the track diagram compressed 
in the ratio of approximately 4 to 1 with 
respect to previous machines. This was 
made possible by removing the push 
buttons and levers from direct physical 
relationship on the panel to the func- 
tions being controlled. A relatively 
small push-button panel is located on the 
desk portion of the machine. By opera- 
tion of the push buttons on this one panel 
the dispatcher may issue controls to any 
desired way station. A logical arrange- 
ment in the numbering of way stations 
and functions permits the operator to 
issue his controls without the need of 
reference to a chart or the memorizing of 
code letters or numbers. 

The material which has been presented 
thus far applies mainly to systems for the 
handling of traffic on single-track lines 
of the Western and Southern parts of our 
country. These same general principles 
have been applied in the reduction from 


Fig. 5. New concept in traffic control machines 
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two tracks to a single track on some rail- 
roads. In some other instances there are 
reductions from three and four tracks to 
two tracks. 

When the broad principles of CTC are 
applied to a railroad with large traffic 
density, such as the rapid transit lines in 
New York City, it is necessary to make 
important changes in the code control and 
indication system. Because of the rela- 
tively short distances involved there is 
not the same urgency to place many 
facilities on one pair of wires. There is 
the requirement for the transmission of 
more information and in less time. This 
is accomplished by using systems with 
greater speed and fewer stations. The 
latest systems permit controls to the field 
simultaneously with indications to the 
office. 

This brief reference to the systems on 
rapid transit lines is made to emphasize 
the range of application of the basic 
principle of directing trains by signal indi- 
cations and to point out that systems 
have been developed to implement this 
principle whether it be on a small stretch 
of track in a remote area, hundres of miles 


of medium traffic density, or on the dense 
traffic of rapid transit lines. 

In the advancement of these develop- 
ments there were always definite goals set 
for the purpose of steering the develop- 
ments along logical lines. One important 
goal that is being taken into account at 
present, and in part guiding development 
programs, is the crewless train. It is ob- 
vious that CTC as discussed in this paper 
constitutes one of the important com- 
ponents in the operation of crewless trains. 
The system transmits the desiresof the dis- 
patcher to the signal system and trans- 
mits information regarding the signals, 
switches, and trains back to the dis- 
patcher. By the techniques used in cab 
signaling, the control information from 
the dispatcher may be applied to the rails 
and transferred inductively to the cab of 
the locomotive. Relays on the locomotive 
controlled in this manner by the dis- 
patcher may have contacts in the loco- 
motive control circuits. Thus the dis- 
patcher would have remote control of the 
locomotive and safety of operation would 
be provided by the way side and cab- 
signaling equipment. 
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may be quite near when crewless tra 


Conclusions 


The material presented in this p; 
and extensive practical experience } 
various installations have led to a x 
ber of conclusions. 


1. The direction of trains from a ce 
control point by control of signal 
tions has advantages over previous 
ods whether the principles be appl: 
light- or heavy-traffic density railro 
cluding rapid transit lines. ; 


2. In general a single-track railroad 
this method of control will have ap 
mately 80% of the traffic capacit 
double-track railroad without this me 
of control. | 


3. After CTC has been installed | 
single-track railroad, the improvem 
the movement of a freight train is i 
by a saving of approximately one minu 
each mile of length of the installation. 


4. The principles of CTC are being 
in some foreign countries to expa 
traffic capacity of existing railroads 
being included in the construction pre 

of new railroads. ; 


5. From a practical standpoint the 


operated on some rapid transit lines. 
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